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Since  nonradiativc  processes  affect  such  critical  parameters  in  solid-state  laser  operations  as  linewidih  of  optical 
trail  .itions  and  energy -level  exchange,  an  understanding  of  the  coupling  between  radiative  and  nonradiativc  proccsse 
is  important.  This  document  reports  upon  one  aspect  of  this  interaction  by  developing  an  internal  phonon  theory 
for  the  thermal  dependence  of  linewidth  in  Nd^^YAG.-^ 

The  approach  to  the  problem  begins  by  first  reviewing  tlie  energy  levels  and  wave  functions  of  rare-earth  tons  in 
crystalline  hosts.  This  esiabli.shcs  notation  and  defines  the  wave  lunctions  which  will  be  used  latci  in  noin.idia 
tive  transition-rate  calculations.  The  allowed  radiative  transitions  arc  then  brietly  considered,  primarily  to  introduce 
radiative  selection  rules  which  will  be  used  later  to  determine  the  allowed  lattice -assisted  vibronic  transitions  (which 
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will  determine  not  only  the  spectral  range  but  also  the  validity  ol'internal-mode  interactions).  The  rigorous  theory  of 
phonons  within  a complicated  crystal  is  then  developed  group-theoretically,  analyzing  the  entire  vibrational  system. 
The  coupling  between  the  lattice  and  the  active  ion  is  subsequently  considered  with  the  single  phonon  process  investi- 
gated rigorously.  A nearest-neighbor  approximation  is  then  presented  and  the  internal-mode  formalism  is  developed. 
Here  the  parameters  which  determine  the  internal  phonon  interactions  with  the  rare-earth  ion  are  explicitly  formu- 
lated and  the  number  of  interacting  modes  predicted.  It  is  the  prediction  of  the  number  of  interacting  modes 
which  will  make  the  internal  mode  approximation  valid  (via  the  vibronic  spectra). 

After  completing  the  formal  derivation  of  the  direct  phonon  transition  rate,  the  two-phonon  Raman  transition 
rate  is  derived,  simplifying  the  expression  where  possible.  The  internal-mode  approximation  is  then  considered  in 
detail  and  the  vibronic  transition  rate  is  derived  to  determine  the  range  of  interacting  phonons  in  the  spectrum.  A 
uniform  phonon  approximation  for  the  “effective”  modes  is  then  proposed,  which  will  determine  not  only  the  range 
of  interacting  modes  but  also  the  dominant  broadening  mechanism.  Before  applying  this  approximation  to  the 
broadening  mechanisms,  the  Debye  linewidth  theory  as  presented  by  Yen.  Scott,  and  Schawlow  is  reviewed.  This  is 
done  primarily  to  place  the  proposed  theory  in  proper  perspective.  Alscr,  this  furnishes  a theoretical  and  historical 
base  to  the  effective-mode  theory  to  be  presented.  Derivation  of  the  equations  which  will  determine  the  functional 
form  of  the  thermal  variation  of  linewidth  is  pursued.  Here  the  effects  of  all  coupled  modes  are  summed  and  analytic 
expressions  are  derived  for  each  of  the  broadening  mechanisms.  Also,  by  employing  the  effective  phonon  spectrum, 
a dominant  broadening  mechanism  can  be  predicted  and  the  linewidth  expressions  greatly  simplified. 

With  this  formalism  developed,  the  comparison  of  the  predicted  linewidth  shape  with  the  experimental  data  is 
investigated  (data  are  available  on  transitions  between  multiplets  of  the  and  41]  and  4lcj^-)).  Only 

normalization  to  room  temperature  (determining  a single  multiplicative  parameter)  and  determination  of  a low- 
temperature  intercept  (dependent  upon  sample  quality)aie  necessary.  The  quality  of  the  fit  to  the  experiment  is 
generally  good  and  within  the  experimental  accuracy.  Tlie  sources  of  error  (especially  at  elevated  temperatures)  are 
discussed  and  are  consistent  with  those  predicted,  ie,  second -order  anharmonic  effects  arc  expected  to  degrade  the 
agreement  at  high  temperatures  since  the  phonon  theory  employed  is  completely  first  order.  Also,  a breakdown  of 
one  of  the  approximations  docs  occur  due  to  the  simplification  made  for  the  effective  phonon  spectra,  wherein 
one  set  of  energy  levels  exhibits  an  extraordinarily  large  direct  phonon  coefficient.  This  is  not  considered  critical  as 
inclusion  of  the  neglected  mechanism  achieved  very  good  agreement.  Possible  limits  to  the  validity  ot  the  approxima- 
tions made  were  indicated,  and  the  trade-off  of  simplicity , predictivity.  and  consistency  is  significant. 
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OBJECTIVE 


Since  nonradiative  processes  affect  such  critical  parameters  in  solid  state  laser  operations 
as  line  width  of  optical  transitions  and  energy  level  exchange,  an  understanding  of  the  coupling 
between  radiative  and  nonradiative  processes  is  important.  This  technical  report  investigates 
one  aspect  of  this  interaction  by  developing  an  internal  phonon  theory  for  the  thermal  depen- 
dence of  linewidth  in  Nd"*"^  YAG. 


RESULTS 

A simple,  compact  phonon  mode  theory  predicting  the  shape  of  linewidth  as  a function 
of  temperature  has  been  developed.  Only  normalization  to  room  temperature  (determining  a 
single  multiplicative  parameter)  and  determination  of  a low-temperature  intercept  (dependent 
upon  sample  quality)  are  necessary.  The  quality  of  the  fit  to  the  experiment  is  generally  good 
and  within  the  experimental  accuracy.  The  sources  of  error  (especially  at  elevated  temperatures) 
are  discussed  and  are  consistent  with  that  predicted,  ie,  second  order  anharmonic  effects  are 
expected  to  degrade  the  agreement  at  high  temperatures  since  the  phonon  theory  employed  is 
completely  first  order.  Also,  a breakdown  of  one  of  the  approximations  does  occur  due  to  the 
simplification  made  for  the  effective  phonon  spectra,  wherein  one  set  of  energy  levels  exhibits 
an  extraordinarily  large  direct  phonon  coefficient.  This  is  not  considered  critical  as  inclusion 
of  the  neglected  mechanism  achieved  very  good  agreement.  Possible  limits  to  the  validity  of  the 
approximations  made  were  indicated,  and  the  trade-off  of  simplicity,  predictivity,  and  consis- 
tency is  significant. 


RECOMMENDATIONS 

Expand  this  work  to  include  other  trivalent  rare  earths  in  other  crystalline  hosts. 
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INTRODUCTION 


The  physics  of  solid-state  laser  materials  is  well  understood.  Both  time-independent 
and  time-dependent  problems,  such  as  line  positions  and  rate  equations,  have  been  investijiated 
thoroughly.  Interactions  of  the  crystal  lattice  with  an  impurity  ion  is  also  fairly  well  under- 
stood, although  not  nearly  as  well  as  atomic  phenomena.  What  is  investigated  much  less  freq- 
uently is  the  quantum  electronics  of  the  nonradiative  processes,  especially  as  it  applies  to  the 
dynamic  operation  of  solid-state  lasers.  Since  nonradiative  processes  affect  such  critical  para- 
meters as  linewidth  and  energy-level  exchange,  an  understanding  of  the  coupling  between  radia- 
tive and  nonradiative  processes  is  very  important.  The  purpose  of  this  report  will  be  to  inves- 
tigate one  aspect  of  this  interaction  by  developing  an  internal  dependence  of  linewidth  of 
trivalent  Neodymium  (Nd'*’'^)  in  Yttrium  Aluminum  Garnet  (YAG). 

The  system  to  be  studied  is  simply  an  example  of  the  general  system  of  a rare  earth  ion 
replacing  a constituent  ion  in  a perfect  crystal  (Nd"*"^  replacing  Y'*’^  in  YACi  I.  What  is  inves- 
tigated about  this  system  is  the  interaction  between  the  impurity  ion  and  the  rest  of  the  lailKe 
and,  in  particular,  how  the  radiative  transitions  are  affected  by  nonradiative  processes.  The 
formulation  of  the  coupling  between  the  active  ion  and  the  surrounding  lattice  will  determine 
how  the  lattice  (ie,  phonons)  affects  the  rare  earth  ion.  It  is  this  interaction  that  broadens 
energy  levels  and  spectrally  broadens  optical  transitions.  When  the  temperature  dependence 
for  the  phonon  occupation  is  included,  the  thermal  dependence  of  linewidth  ma\  be  derived 

In  order  to  find  the  dominant  broadening  mechanism  as  well  as  the  functional  lorm  ol 
these  interactions,  the  spectrum  of  those  phonons  affecting  the  active  ion  must  be  vietermined 
In  previous  treatments,  best  illustrated  by  Yen,  Scott  and  Schawlow'-^"’-'’.  the  Debye  approxi- 
mation was  used,  which  considers  only  acoustic  modes.  This  "external"  mode  theory  i>.  a pooi 
approximation  to  what  actually  occurs,  since  lattice  vibrations  within  the  unit  cells  are  ignored 
The  theory  presented  here  departs  from  the  previous  treatments  by  considering  only  those 
optical  phonons  which  couple  strongly  to  the  impurity  ion;  ie.  those  internal  phonons  result- 
ing from  nearest-neighbor  interactions.  This  leads  to  a theory  much  different  in  intent  as  well 
as  in  functional  form.  The  objective  of  the  Debye-type  theories  is  not  to  predict  the  shape  id 
the  linewidth-vs-temperature  cui-ve.  nor  even  to  formulate  analytic  expressions  for  the  broad- 
ening mechanisms  vs  temperature,  but  rather  to  check  the  validity  of  the  proposed  broadening 
mechanisms  themselves.  The  intent  was  never  the  formul.itii'ii  of  a simple  theory  able  to 
predict  the  shape  of  the  linevvidth-vs-temperatiire  curves  in  trivalent,  rare-earth  dopeil  cry  sUiK. 

This  is  the  intent  of  this  technical  report:  establish  the  internal  mode  formalism  necess.iry  to 
model  the  interactions;  determine  an  approximation  for  the  "effective"  phonon  spectra; 
using  this  approximation,  derive  coin|xict  expressions  for  the  broatleniiig  mechanisms  that 
do  not  reviiiire  numerical  integration  nor  the  determination  of  nonmultiplicative  parameters; 
employ  approximatii  ns  to  simplify  the  rornniliMn  ,md  tletermine  the  dominant  broadening 
mechanism;  aiul,  fiiKilly.  to  attempt  to  prcilicl  the  shape  of  the  linevvidth-vs-temperature 
curves  for  the  transitions  on  which  data  are  available.  The  comparison  to  experimental  results 
will  be  found  to  be  within  experimental  .iccuracy  of  the  ilata.  and  sources  of  error  are  to  be 
examined. 

rite  appro;ich  to  the  problem  begins  by  first  reviewing  the  energy  levels  and  wave 
functions  of  r.ire-earih  ions  m cry  stalline  hosts.  | his  establishes  noi.ilii'n  and  del  ines  the 

' \V  Siiiiili  .iikl  I’  Sorokin , I he  I ,iser . Meduivv  -I  I ill  book  (onip.inv  , t'Hiu 
‘\V\t  Yen . yv  (' Seoii . anil  .\( ' Sili.iw  low  . (’hy  sieal  Review.  I .hi.  .A’T  I . pio4 
' 1)1  Me(  limber  anil  Ml)  Si  uu:o . .loin M.il  ot  Applied  l*liysics.  l(>S2. 
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wave  riinctions  wliieli  will  be  used  later  in  nonradiative  transition-rate  ealeulations.  The 
allowed  radiative  transitions  are  then  brielly  eonsidered.  primarily  to  introduce  radiative 
selection  rules  which  will  be  used  later  to  determine  the  allowed  lattice-assisted  vihronic  trans- 
itions (which  will  determine  not  only  the  spectral  range  but  also  the  validity  of  internal-mode 
interactions).  The  rigorous  theory  of  phonons  within  a complicated  crystal  is  then  developed 
group-theoretically . analyzing  the  entire  vibrational  system.  I'he  coupling  betw'cen  the  lattice 
and  the  active  ion  is  subsequently  considered  with  the  single-phonon  process  investigated 
rigorously.  A nearest-neighbor  appro.ximation  is  tlien  presented  and  the  internal-mode 
formalism  is  developed.  Here  the  parameters  which  determine  the  internal  phonon  interac- 
tions with  the  rare-earth  ion  are  explicitly  forimdated  and  the  number  of  interacting  modes 
is  predicted.  It  is  the  prediction  of  the  number  of  interacting  modes  which  will  make  the 
internal-mode  approximation  valid  (via  the  vihronic  spectra). 

.\fter  completing  the  formal  derivation  of  the  direct  phonon  transition  rate,  the  two- 
phonon  Raman  transition  rate  is  deriveil  and  the  expression  is  simplitied  where  possible.  The 
internal  mode  apjsrc.ximation  is  then  consiilered  in  detail  and  the  vihronic  transition  rate  is 
derived  to  determine  the  range  of  interacting  phonons  in  the  spectrum.  A uniform  phonon 
approximation  for  the  ''etfective'  nuules  is  then  proposeil  which  will  determine  not  only 
the  range  of  interacting  modes  but  alsci  the  dominant  broadening  mechanism.  Before  appK  - 
mg  this  appro.ximation  to  the  broadening  mechanisms,  the  Debye  linewidth  theory  as  pre- 
sented by  Yen,  et  al.  is  reviewetl.’  This  is  done  prim.iriK  to  place  the  proposed  theor\'  in 
proper  perspective.  Also,  this  furnishes  a theoretical  ,md  historical  base  for  the  effective-mode 
theory  to  be  presenteil.  Derivation  ol  the  ei|uations  which  will  determine  the  functional  lorm 
of  the  thermal  variation  of  linewidth  is  now  pursuetl.  Here  the  effects  ol  all  coupled  modes 
are  summed  and  analytic  expressions  are  derived  for  each  of  the  broadening  mechanisms.  .Also, 
by  employing  the  effective  phonon  spectrum,  a dominant  broadening  mechanism  can  now  be 
predicted  and  the  linewidth  expressions  greatly  simplified.  W ith  the  formalism  developed,  the 
comparison  of  the  predicted  linew  idth  shape  w ith  the  experimental  data  is  investigated.  Only 
normalization  to  room  temperature  (determining  a single  multiplicative  parameter)  and 
rletermination  of  a low  temperature  intercept  (dependent  on  sample  (luality ) are  necessarx  . 

I he  (|uality  of  the  fit  to  the  experiment  is  generally  good  and  within  the  experimental 
accuracy.  The  sources  of  error  (especially  at  elevateil  temperatures)  are  discusserl  am.1  are 
consistent  with  those  preilicteil.  ie.  second-order  anharmonic  effects  are  expected  to 
degraile  the  agreement  at  high  temperatures  since  the  phonon  theory  employ  ed  is  com- 
pletely f irst  order.  Also,  a breakdown  of  one  of  the  approximations  does  occur  tiue  to  the 
simplification  made  for  the  effective  phonon  spectra  wherein  one  set  of  energy  levels 
exhibits  an  extraordinarily  large  direct  phonon  coefficient.  T his  is  not  considered  critical 
as  inclusion  of  the  neglecteil  mechanism  achieved  very  good  agreement.  Possible  limits  to 
the  validity  to  the  approximations  made  were  iiulicated.  ,md  the  trade-off  of  simplicity . 
predictivity  . and  consistency  is  significtmt . 


SI  MM  ARV  Ol  BRO  ADIM\(.  MK  H \MSMS 

If  one  IS  to  formulate  the  effects  cif  the  phonon  interactions  .ind  derive  the  phomm- 
photon  equations,  one  must  first  understand  non-raili.itive  inter.ictions.  Ihis  will  lead  It'  an 
undersi.mding  ol  the  thermal  depeiulence  ol  linewidth.  l et  us  now  brielly  review  some  of 
the  physics  involved.  The  width  of  a transition  (linewidth)  is  affected  by  four  mechanisms, 
each  having  a different  lemper.iture  depeiulence.*  first  is  the  mech.inisin  (d  random  micro- 
sco(nc  strain  within  the  host.  Since  the  statistic.il  varialK'ii  in  .1  cryst.il  lield  im  the  impurity 


(1 


ion  results  in  a variation  of  tlie  Stark  splitting,  we  expect  a Ganssian  broadened  line  that 
should  be  temperature  independent.’’'’ 

The  widths  of  radiative  transitions  (although  ions  vibrate  rapidly)  exhibit  no  Doppler 
broadening  because  the  amplitude  of  vibration  is  so  small  If  ions  have  very  short  lifetimes  in 
a particular  state,  then  Heisenberg  broadening  will  give  rise  to  homogeneous  Lorentzian  lines 
having  widths  proportional  to  the  total  transition  rates  of  each  level.  Consider  the  typical 
energy  level  scheme  encountered  with  a rare-earth  ion  in  a crv'stal  (figure  I ).  Here  we  find 
groups  of  crystal  field  split  Stark  levels  for  each  spin-orbit  multiplet  which  are  characleri/ed 
by  2S-M  [ [ (energy  levels  in  rare  earths  will  be  examined  later).  Consider  the  state  The 
total  transition  rate.  . both  radiative  and  uonradiative.  upward  or  downward  to  Stark 
levels  within  the  J manifohl  or  to  other  J manifolds,  determines  the  Inmiogeneous  lincwiillh 
( full  width,  half  maximum  ) 


Since  the  radiative  lifetimes  are  typically  on  the  order  of  l()~-  to  10“^’  second,  we  can  be 
sure  that  this  is  negligible  with  regard  to  the  observed  linewidths.  An  understanding  of  the 
mechanisms  of  uonradiative  energy  level  interactions  will  thus  lead  to  the  thermal  variation 
of  linewidth  and  to  an  understanding  of  the  energy  level  exchange  iluring  the  operation  of 
a laser. 


Al  RADIATIVE  TRANSITIOM 
J' 

Bl  SINGLE  PHONON  TRANSI 
TION  , 

Cl  MULTIPHONON  TRANSinON 


2S+1 1 
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If  there  are  phonons  in  the  lattice  resonant  with  the  spacing  of  two  ionic  energy- 
levels,  then  the  lattice  may  couple  with  the  ionic  states  and  a direct,  single  phonon,  non- 
radiation transition  may  occur  either  by  phonon  absorption  or  emission  (figure  2).  We  will 
show  later  that  and  (the  downward  and  upward  direct  transition  rates)  behave  like: 


[n(o;)+  HD'(cu) 
oc  n(cj)  D'  (co), 

where  n(u>)  is  the  occupation  number  for  phonons  of 
number  is 


(1) 


frequency  oj,  the  thermal  occupation  j 

i 


n(cu)  = [exp(hco/kT)  - 1 ] *, 

and  D(w)  is  the  effective  density  of  states  for  the  ion-lattice  system.  The  energy  level  spacing 
of  the  system  of  Stark  levels  is  typically  1 0 ' to  1 0^  cm"'  which  results  in  transition  rates 


DIRECT  UPWARD  RAMAN  UPWARD  MULTIPHONON  TRANSITION 

TRANSITION  TRANSITION  (P  PHONONS  OF  ENERGY  A) 


• virtual  or  intermediate  STATES 

figure  2.  Illusiratioiis  ol  the  phonon  broadening  mechanisms  involving  nonradiative  transitions. 
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from  l v'i>  of  about  10~*  ' second  wliidi  results  in  levels  ~3  cnr*  wide.  1 lieretore.  it  Ilie 
transition  rates  are  summed  tor  interactions  witli  all  the  other  levels,  the  total  iransiiion  rate 
from  |i//j>  can  be  determined  and.  hence,  the  width  of  If  the  same  is  done  tor  vq.-'*- 

then  the  width  of  the  optical  transition  from  |i^j>  to  niay  be  determined  I lie  plionon 

depeiulence,  and  tlierefore  the  temperature  dependence,  ot  that  hnewidth  is  then  knovui 

If  the  level  spacing  is  such  that  hco  > where  cC|ii;,x  i''  >l'e  laiecst  irec|uenc\ 

phonon  supported  by  the  lattice,  then  a higher-order  process  must  be  iiuokcd.  one  which 
involves  the  emission  and  absorption  of  many  phonons.  In  general,  tlie  multiphonon  relax- 
ation rate  goes  down  as  the  order  of  the  process,  so  in  the  energy  conser\ation  eipiation 
A X p = li,_,/h.  we  try  to  minimize  p and  still  have  A contained  in  the  phonon  spectrum 
where  p is  the  number  of  phonons  (order  of  the  process),  is  the  energ\  gap.  and  A is  the 
phonon  freciuency  (figure  2).  Thus,  we  would  expect  a simple  expression  of  the  form 

|(n(A)-r  l)l)(A)|P 
cl 

C) 

ln(A)l)  (A)|P  (.3.  0,7.8) 

for  the  downward  and  upward  multi  phonon  transition  rates.  Tlie  probability  of  a more  com- 
plicated process  reciuiring  phonons  of  many  freiiuencies  could  also  occur  because  of  the 
shape  of  t)  (co),  selection  rules,  and  so  on.  Since  W'”P  is  a function  of  the  occupation  number 
n for  phonons  A,  it  lias  a temperature  dependence  ami  is  in  contradiction  to  the  article  by 
Yen.-  Approximate  temperature  independence  is  achieved  only  when  n is  small  compared 
to  1 in  eiiuation  (2).  Because  A is  large,  n is  small  at  all  but  the  highest  temperatures.  It  is 
important  to  note  that  it  is  this  process  whicli  empties  the  thermal  level  to  the  ground 
state  in  Vi’*'-''  'i'.ACi.  The  emission  of  two  '-600-cm~'  phonons  account  for  the  ().,x-/as  deca\ 
time  from  the  thermal  state,  riicrefore.  while  the  tmiltiphonon  process  contributes  little  to 
the  linewidtii  broailening,  it  is  very  important  in  the  exchange  of  energy  levels  (especialh 
the  thermal  level  deiiopulation ). 

1 1 remaining  broadening  mechanism  involves  a second-order  process  which  pre- 
dominates at  higher  temperatures,  the  Kamaii-lik  .•  'Cattei  u g of  phonons  from  impurii  v ani' 
(eg.  a Nd^  ' ion  replacing  the  ion  in  Y<\(ii.  1 his  can  be  thought  of  as  an  absorption  oi 
one  phonon  and  the  emission  of  another  with  the  differenee  frei|uenc\  correspoiuling  to 
the  energ\'  gap  of  the  transition.-  I'or  a deca\  |irocess  we  find  the  transition  rate  between 
levels  separated  by  Al  (cos  - c<j  | = Al  h ) is  just 

\V*^  a D'tcc  I )n(co  I I X D'lcc  s )|  nice  1 ) + I I . ( I 

where  (Ct  and  cjs  are  the  frequencies  ot  the  phoiunis  ( t ieure  2).  We  can  see  that,  while  at 
low  temperatures  n - 0 and  so  iloes  W . at  high-temperatures  W'  can  become  the  predom- 
inant efiect.  the  n's  increase  md  many  more  phonons  can  contribute  since  ill  th.it  i'  requir  'd 
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is  that  C0-)  - W|  = ^L'h.  Therefore,  tlie  total  transition  rate  is  summed  for  all  interactinj: 
phonons  in  the  spectrum. 

liven  from  this  brief  description,  the  importance  of  knowing  the  detailed  density  of 
states,  D (cu),  is  obvious.  For  e.xample.  we  know  that  the  phonons  wliich  cause  Raman-like 
rela.xations  are  those  where  the  difference  between  two  peaks  in  D (oj)  corresponds  to  the 
transition  energy.  Dtco)  is  also  important  in  determining  the  order  of  the  multiphonoti  pro- 
cess, and  so  on.  The  use  of  the  Debye  acoustic-phonon  model  is  common  in  deriving  line- 
width  expressions  but  is  actually  a very  poor  approximation  to  what  occurs.  In  the  follow- 
ing sections  a more  correct  expression  will  be  used  which  limits  the  phonon  spectrum  to 
onl\  those  optical  phonons  which  couple  strongly  to  the  impurity  ion.  This  results  in  a 
reduction  in  the  importance  of  the  direct  process  and  simplifies  the  Raman  expression. 
Finally,  a linewidth  expression  is  derived  which  fits  the  experimental  results  to  within  a few 
per  cent,  within  the  stated  experimental  accuracy. 


TRIVALENT  RARE  EARTHS  IN  CRYSTALLINE  HOSTS 
ENERGY  LEVELS  AND  WAVE  EUNCTIONS 


The  treatment  of  a rare-earth  impurity  ion  located  in  a crystalline  host  is  very  com- 
plicated. In  order  to  begin  to  understand  the  quantum  electronics  of  selection  rules,  phonon 
transitions,  and  ion-lattice  effects  in  general,  a good  understanding  of  ionic-energy  levels 
and  wave  functions  is  required.  This  section  reviews  the  necessary  physics  and  establishes 
notation  to  be  used  later.  The  system  of  a Nd”^^  ion  located  in  a Y'*'^  site  in  YAG  is  used  as 
an  example;  the  formalism  necessary  to  handle  general  ion-host  combinations  is  similar. 

Atomic  Nd  is  number  60  in  the  periodic  table  and  has  a configuration  of  XENON 
-t-  (df)'^  (.Sd)^  (6S)-.  The  normal  oxidation  state  of  rare  earths  is  the  trivalent  state  which 
results  in  a ground-state  configuration  for  Ny^^  of  XENON  + (4f)^  (5d)0.  This  is  not  the 
only  possible  configuration  but  merely  the  lowest -energy  one.  Other  possible  configurations 
are  (40-  (.^id  ) ' , (4f)-  (6Sl'.  both  having  even  parity.  The  importance  of  higher-energy 
configurations  becomes  clearer  if  parity-forbidden  electric-dipole  transitions  are  investigated. 

Yttrium  Aluminum  Garnet  has  a chemical  composition  Y^  Ah  (.\1  Oqt-^  which  is 
composed  of  body-centered  cubic  cells.  The  arrangement  of  the  80  ions  which  compose  the 
unit  cell  is  described  by  the  space  group  symmetry.  The  local  field  (or  that  which  is  seen 
by  the  Ntl'*'-^  ion)  is  described  by  the  point  group  symmetry  ; in  this  case.  D-i  (or  222  in 
international  notation).  D-i  symmetry  is  that  of  three  mutually  perpendicular  rotations  of 
rr  (note  figure  3).  Notice  that  the  crystal  does  not  have  inversion  symmetry  and,  therefore, 
parity  will  not  he  a good  c|uantum  number.  When  this  ion-host  combination  is  used  as  a 
laser  material,  the  Nd'’’-''  ions  are  substituted  into  Y'*'-^  sites  at  approximately  one  atomic 
per  ce  lt,  l.arger  doping  levels  reduce  the  validity  of  a single-impurity  ion  approximation  in 
the  lattice  (physically  we  get  ion-ion  effects  which  have  macroscopic  effects  such  as  reducing 
the  radiative  quantum  efficiency  and  inducing  crystal  strain). 

I he  Hamiltonian  of  the  Nd'*'-^  ion  in  the  crystalline  site  may  be  written  as 
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Figure  3.  Character  table  of  the  irreducible  representations  for  the  crystal  double  group 
Di  and  the  stereographic  projection  of  the  Ds  point  group. 
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The  first  term  is  the  kinetic  energy  of  tlie  N electrons,  the  seeoni.1  the  electron-nueletir 
attraction,  the  third  tlie  electron-electron  cotilonih  repulsion,  the  fourth  the  spiti-orhit 
interaction,  and  the  hist  the  static  effect  of  the  lattice.  Here,  we  also  have  F|^  as  the 
position  of  the  k'^’  electron,  f(r|^)  as  the  spin-orbit  coefficient,  aiul  N = ,37.  7.  = 60.  Since 
direct  determination  of  the  eigenvalues  of  II  i-.  impossihli. , > .^eessive  applications  of  firsl- 
order  degenerate  perturbation  theory  is  uscil.  In  verv  simple  terms,  we  divide  up  the  ionic- 
parts  of  1 1 as  follows: 


1 -^^k"^'"-k* 


N N ■>  N 

lies  = - I Vtr^t  + ^ ^ 

k = 1 i>j=l  " I k = I 


"so  " ^ '■k  ■ 'k 
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(in  the  central  field  approximation)  determines  the  electronic  configuration,  (4f)^.  and 
the  starting  wave  functions.  This  configuration  is  364-fold  degenerate  with  respect  to  m^i  and 
m^.  The  electrostatic  perturbation  has  the  effect  of  adding  the  individual  6’s  and  s’s  and 
results  in  wave  functions  described  as  I?  L S Mj^  Mg  > with  L = 

N _ _ N 

^ S = and  ? designating  all  the  other  quantum  numbers  of  the 

k=  I k=  1 

configuration  such  as  the  n and  ^ designation  of  (4f).  Formally,  the  procedure  is  to  find  a 
vector  space  of  functions  composed  of  linear  combinations  of  the  old  symmetrized  wave 

functions  In^  nig  m^  > |m^.,  nig  in^  > • • ■ which  diagonalize  Thus. 
a a b b 

splits  the  364  degenerate  levels  into  different  levels  depending  upon  L and  S (the  terms). 

Facli  term  is  still  (2L  + 1)  (2S  + 1)  fold-degenerate  with  respect  to  Mj^  and  Mg. 

The  application  of  the  spin-orbit  perturbation  results  in  the  vector  addition  of 
1.  and  S ; where  formally  we  say  that  the  terms  are  diagonalized  in  the  space  L S J Mj  >, 

N 

with  J = N (t|^  + S|^).  Now  the  wave  functions  are  only  (2J  + 1 ) degenerate  with 
k=  I 

respect  to  Mj. 

The  validity  of  applying  the  perturbations  in  this  order  (the  LS  coupling  scheme)  is 
confirmed  spectroscopically  ( figure  4).  where  we  find  that  ~1 0,000  cm”' , ~ 1 .000 

cm~' . and  H^.j.y^j.,1  ~ 100  cm“' . Thus,  we  expect  a given  term,  labelled  by  L.  to  be  split 
by  spin-orbit  coupling  into  a close  group  of  levels  which  are  still  characterized  by  L and  S (to 
a good  approximation)  and  distinguished  by  different  values  of  J,  noted  by  2S+1  j_j  Ji^e 

reduction  of  the  (2J  + 1 ) degeneracy  is  accomplished  by  the  static  Stark-effect  term.  I^^•rv;;lal• 
The  fact  that  the  ion  interacts  only  weakly  with  the  crystal  lattice  is  due  to  the  small  spatial 
extension  of  the  (4f)  wave  functions.  Likewise,  this  will  be  responsible  for  narrow  spectra 
from  transitions  lying  wholly  within  the  (4f)^  configuration  (whereas  interconfigurational 
(4f)-^  ^ (4f)-  5d  transitions  are  characteristically  much  broader). 

I he  application  of  the  H^TVstal  Pe'fturbation  term  turns  out  to  be  much  more  diffi- 
cult. .Mthough  a brute-force  diagonalization  within  the  (2J  + 1 ) ionic  manifold  is  certainly 
possible,  mathematical  inconveniences  render  it  impractical.  Instead,  by  using  the  essential 
property  of  the  crystal,  namely  symmetrv’,  and  the  group  theory,  the  determination  of  the 
(|ualitative  splitting  due  to  the  crystal  potential  is  greatly  simplified.  More  important  though, 
the  transformation  properties  of  the  wave  functions  are  detemfined  which  allow  phonon 
selection  rules,  ion-lattice  effects,  and  the  like  to  he  systematically  studied. 

In  't'ACi.  with  the  Nd'*’-^  ion  located  in  the  low  site  symmetry  Dt,  all  except  the 
residual  two-lold  Kramer  degeneracy  is  removed.  The  character  table  of  the  crystal  double 
group  derived  for  the  nonintegral  .'’s  is  given  by  figure  3.'* 

Since  the  purel\  ionic  II  is  spherically  symmetric,  we  find  that  the  symmetiy  of  the 
total  II  is  that  of  thisgroup  theoretical  notation,  it  is  determined  ' .lat  all  wave 

functions  are  described  by  a two-dimensional  F representation  within  the  crv  :al  double 
group  l)s.  For  example,  the  main  laser  line  is  the  transition  (4f-^/s)->  F (41|  j ^s);^  F'  (the 
2 means  the  secotui  Stark  level  from  the  bottom).  It  is  only  because  of  the  low-point  group 
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Kigurc  4.  Observed  energy  levels  of  the  trivalent  rare-earth  ions.' 

s\  miiK  try  th;it  vve  have  all  energy  levels  transforming  the  same;  in  general,  different  wave 
futietiotis  corresponding  to  different  states  will  transform  tlifferently  under  operations  of 
the  group  and  retain  degeneracies. 

.Although  grotip  theory  is  very  powerftil.  it  cannot  give  the  amount  of  Stark  split- 
ting (matrix  elements  must  still  be  lotmd).  Ihe  techniciues  for  calculatitig ciuantitative 
splitting  are  not  discussed  here.'’  The  energy  levels  of  Nil'*’-^  \ .ACi  are  given  in  figures  5 and 
b. 


ALLOVVFI)  RAOIATlVt  TR  ANSITIONS 

Transitions  between  these  levels  are  responsible  for  the  operation  of  these  maleriaK 
as  lasc.  >.  While  the  calculation  of  radiative  transition  rates  for  rare  earllis  is  ;ictively  puisued.'"  •" 


only  a brief  review  will  be  presented  here.  The  observed  spectra  of  RE  ions  consist  mainly  of 
electric  dipole  transitions’’'®  which  are  governed  by  matrix  elements  of  the  type  <f  I - e F | i>. 

If  the  initial  and  final  states  I i>  and  I f>  have  the  same  parity  (as  with  crystals  having  inversion 
symmetry  and  for  transitions  taking  place  wholly  within  the  4f  shell)  the  electric  dipole  transition 
is  rigorously  forbidden.  Only  if  there  are  odd  parity  terms  in  the  static  crystal  potential  or  if  there 
are  odd  pari^  vibrations  to  admix  opposite  parity  states  (eg,  wave  functions  derived  from  the 
excited  (4f)-5d  configuration)  can  electric  dipole  transitions  occur. We  will  assume  here 
that  only  the  odd  terms  of  the  crystalline  potential  can  cause  mixing  of  parity  states  and  that 
the  vibrational  contribution  is  small  (later  we  will  use  the  weak  vibronic  emission  and  absorp- 
tion spectra  to  determine  the  range  of  interacting  phonons). 

As  stated  before,  the  total  electronic  Hamiltonian  is  H^i  = H(Yeg  ion 

M^.  = -e  Sj  V,(ri.0i,0i), 


V^,  (rj.  Oj,  0j)  being  the  crystal  field  potential  at  the  i-th  (4f)  electron  (in  our  example  the  sum- 
mation extends  over  the  three  (4f)  electrons  of  Nd"*"^).  It  is  often  convenient  to  expand  H^. 
in  spherical  harmonics 


'•c=  :i 

n m = 
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(6) 


where  the  coefficients  are  constants  determined  by  the  position  of  the  neighboring  ions 
around  the  central  RE  ion  (ie,  point-group  symmetry),  r'|Vefers  to  the  radius  of  the  i-th  (4f) 
e]ectron  and  the  (ffj.  0,)arc  spJicrica)  harmonics.  This  potential  may  be  separated  into 
even  and  odd  parity  parts  (with  n even  and  odd  respectively)  and  written  as 

11^,  = H^.  (even)  + 11^.  (odd). 


In  the  calculation  of  line  positions,  only  the  even  terms  of  M^.  are  important  and.  therefore, 
the  othl  spherical  harmonics  are  ignored.  Since  the  Stark  cplitting  involves  matrix  elements  such 
as  < VtVee  ion  I ^ *^free  ion  consider  terms  in  11^,  with  n < 6 (since  the  1 

electrons  with  (’  = 3 cannot  connect  states  with  n > 3 + 3 = b).'"' 

Analogous  results  are  obtained  when  calculating  transition  rates.  Here  we  w ill  reipiire 
matrix  elements  between  admixeil  states  such  as  < f | - e r | x >.  with 

<d  I H . (odd)  I i > l)3> 

lx>=li>+^^ . (7) 

' i - 'd 

where  the  summation  extends  over  all  excited  opposite  parity  states  (3-  These  states  are 
derived  from  higher  energy  configurations  than  4f".  such  as  4f'’’~'  5d  (note  figure  "i.  Since 
11^,  is  limited  to  terms  with  n < (n  this  imposes  a selection  rule  which  tends  to  decrease  the 
line  intensity  with  increasing  .\J  up  to  AJ  = |h,>,  is  one  reason  wh\  the  41  , ^4I|<;^-, 

transition  (AJ  = (i)  is  approximalely  two  orders  ol  magnitude  less  intense  than  "^EA  s ^ “^l  I I zs.  ~ 


'^JD  Judd.  Physical  Rcucw.  i:7.7(d),  PU,: 

"H(  Wyhiiiiiiic.  joum.il  i)l  riuMiiisUs  .uid  l’h>  siss,  72 , lUuO 
l,\  Kuiimiisu-iii  .Mid  It  (iciiM.  IMiysic.il  Review.  I.!(i.,\7||.  l'l(,4 


l.x' 


4 


l93l 

ImJ 


Ce  Pr  Nd  Pm  Sm  Eu  Gd  Tb  Dy  Ho  Er  Tm  Yb  Lu 


(igiiro?,  lTicrg>  IcvL-Uliagi  am  ol  llio  4l"  ( w hiu- ) and  41 (t’lack  I 
toiitigiiralioii  i>t  irivalonl  ran.' -earth  ion.-. 


Here  apain.  while  tlie  perUirhalive  Ireatmeiil  ei\es  i|u.inlitatlve  results,  the  croup  Iheor 
etica!  approach  is  useii  to  iletermme  simple  selection  rules  (this  methoil  only  fituis  torbidden 
transitions,  not  allowed  onesl.  .Since  the  ilirect  product  ol  I'  ( i >)  X I"  1 1 r>)  when  reduced 
in  the  site  symmetry  l)s  contains  irreducible  represent. itions  ol  the  tr.insition  moment  r (odd 
ptirity  crystal  terms  arisinc  from  V^-.  spherie.il  li.irmonKsi.  no  electric  dipole  transi 

tions  are  forbidden  lie.  no  distinction.  ;it  letist  Irom  croup  theor\  , .is  to  tt  ;ind  o electric  dipole 


radiation).'^’**  It  is  again  only  the  very  low  site  symmetry  of  D->  that  prevents  selection  rules; 
in  general,  different  selection  rules  are  obtained  for  transitions  between  states  characterized  by 
irreducible  representations. 


SINGLE  PHONON  INTERACTIONS 


PHONONS  IN  YAG 


Consider  the  system  of  a single  rare-earth  (RE)  ion  replacing  an  ion  in  a perfect  cr>  stal. 
The  Hamiltonian  for  this  system  is 

H = -H  Hp (8) 

where  is  the  electronic  Hamiltonian  of  the  RE  ion  in  the  perturbing  crystal  field  of  the 
other  ions  situated  at  their  equilibrium  positions: 

Hei  - jQP  + Hg, 

Hg  is  the  crystal  potential  energy,  Up  is  the  phonon  Hamiltonian,  and  Hgp  describes  the  inter- 
action between  the  electronic  state  and  the  lattice  phonons.  The  zero-order  wave  functions 
I i//  >,  which  are  solutions  of  the  equation 

Hg,  I i//>  = E I !//>, 

are  characterized  by  the  frte-ion  quantum  numbers  fLSJ,  an  irreducible  representation  of  the 
RE  ion  site  symmetry  group  G^  (ie,  the  point  group  symmetry  in  the  vicinity  of  the  RE  ion), 
and  a row  of  this  representation  (to  specify  each  wave  function  within  a degenerate  represen- 
tation). 

Consider  now  the  second  two  terms  in  H contained  in  equation  (8).  The  phonon 
Hamiltonian  Hp  is  invariant  under  the  operation  of  the  entire  space  group  G.  The  vibrational 
wave  functions  are  characterized  by  irreducible  representations  of  this  group.  If  a crystal  has 
N unit  cells  and  n ions  per  unit  cell,  we  expect  3Nn  possible  modes  of  vibration.  The  size  of 
the  matrix  to  be  diagonalized  is  reduced  from  a 3Nn  X 3Nn  matrix  to  a 3n  X 3n  matrix  b\ 
considering  translational  symmetry;  all  possible  values  of  k being  contained  within  a Brilloum 
zone  (BZ).  There  now  remain  3n  branches  of  the  phonon  spectrum  to  be  obtained  with  each 
branch  being  defined  by  an  irreducible  representation  of  G.  By  convention,  the  3n  - 3 iikhIcs 
which  ue  not  have  zero  frequency  at  k = 0 are  called  optical  modes  and  those  with  ccik  = Oi  = 0 
are  acoustical  modes.  It  should  be  noted  that  the  degeneracy  of  particular  branches  ma\  be 
lifted  along  directions  of  lower  symmetry  of  vector  k than  k = 0.  In  the  case  of  ^ .\1  s ( ,M 

the  space  group  is  la  3d  (0|^)*^,  a body  centered  cubic  Bravas  lattice.  There  are  80  transla- 
tionally  invariant  sites  (n  = 80)  in  YAG,  and  by  reducing  the  reducible  representation  of  ionic 


'' J(  I’i.iIIkm  . ,\lomie  bnergy  Levels  in  Crystals,  National  Bureau  of  Standards  Monograph  l‘>.  Fehniats  I'tol , 
US  (ioveinment  I’miting  Office.  Washington.  DC 
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displacements  in  each  of  the  four  types  of  sites  ( Y transforms  like  Dt,  O like  C ] and  the  two 
types  of  A1  sites  like  S4  and  S^),  the  symmetry  of  all  possible  phonons  at  «=  0 can  be  found 
(figure  8): 


5A]y  + 3A|j,  + 5Asu  + 5Asg  + lOE^j  + 8Eg  + 14T  jg 
+ 18T  ]u  + I'fTig  + 16Ts^j. 


(9) 
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figure  8.  The  character  table  of  the  reducible  representations 
generated  by  arbitrary  displacements  5 within  YAG. 

The  above  constitutes  a total  of  240  modes  (ie.  3 X 80  branches  of  the  cu(k)  dispersion  curve). 

By  assuming  the  harmonic  appro.ximation.  the  phonon  Hamiltonian  Hp  becomes  the 
sum  of  the  harmonic-oscillator  Hamiltonians  for  each  normal  coordinate; 


V 

k>r 


J.  r2r--  0.-L  - nr  n*r 


<^k 


k7  k7  k7 


(10) 


rite  normal  coordinate  is  characterized  by  a wave  k,  a small  representation  7,  and  a row 
r of  this  representation  (eg.  to  account  for  all  polarizations  within  a T|g,  three-fold  degenerate, 
vibration ).  The  Schroedinger  equation  for  a particular  normal  coordinate  Qi.,^  is 


It"'  1 *>  r *r 

h-  V^r  0),;^  0^;^ 


I <^k7  "k7 


> 


(III 


= 'nk7+3)  '’r^k7'^k7  "k7  >• 
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where  and  contain  no  subscript  r since  different  rows  within  a representation  indicate 
degeneracies.  Therefore,  the  wave  function  of  the  whole  system  of  phonons  is  the  product 


n I Ok;  Hk;  > 

k7r 


(12) 


and  the  corresponding  energy  is 


k7r 


<"k7  +-)tiWkv 


k7' 


(13) 


The  three  components  of  the  displacement  Xj(n;a),  i = 1 , 2,  3,  of  the  ion  in  the  n^^  unit 


cell  are  related  to  the  normal  coordinates  by  the  linear  transformation 


Xj(n;a)  = ^ ( N inQ,)”'  - ( k ; a)  exp  (i  k • R^) Qk^  (12) 

k7r 


where  is  the  position  of  tlie  n**'  unit  cell  and  (k;a)  is  a measure  of  how  much  a normal 
mode  Ok^  affects  the  particular  displacement  Xj(n;a).  In  terms  of  boson  operators. 


r _ nr  / + \ 

■ Y2cuk^  I 


There  are  two  orthonormality  relations  for  the  Cj.y|.'s  given  by 


^ ( i7r  < «)  f'i7'r'  *k.a) 
\a 


V ♦ - * , 

^ ( i-yr  * k • Q • f j-yr  ( )>  • cr  • ~ i^jj  ^cva'- 

7r 


l.et  \’  be  the  total  interaction  Hamiltonian  between  the  RE  ion  and  the  rest  ot  the 
ions  in  the  crystal.  It  can  he  expanded  in  terms  of  the  displacements  Xj(n;a): 


Vo"  ^ 


noi  = 1,2,3 


Xjln.a)  7j  (n;a)  V’  | q + . 


(14) 


The  first  term  Vq  is  the  electrostatic  (Stark)  term  which  is  incorporated  in  H^.|  as  11^..  In  the 
second  term,  the  derivatives  are  evaluated  at  the  equilibrium  positions  of  all  the  ions,  with  the 
Rli  ion  considered  to  be  at  the  origin.  By  using  equation  ( 1 2)  we  can  write  the  interaction 
Hamiltonian  between  the  RE  ion  and  the  lattice  vibration  as 
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ep  = X 
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exp  (ik  • R,^)  Vj  (u;q;I  V 


k7r 

• Oi.I.  = vi  f 


(N  iiIq,)"  - Cj.y^  (k;a) 


Qk;=  vv>kV 

k7r 


(151 


4 


1 


t 

The  potential  gradient  terms  defined  above  operate  only  on  the  RE-ion  electronic  state 

while  0]^^  (containing  a and  a"*"  terms)  operate  on  phonon  states  only.  This  separation  in  V^p 
is  valid  under  the  assumptions  stated  before,  namely:  V^p  « and  Vg  « ie,  the 
mixing  of  electronic  and  vibrational  states  is  small.  Since  we  have  assumed  only  a single  RE  ion 
interacting  with  the  lattice  phonons,  this  implies  that  Vg  is  not  invariant  under  the  full  space- 
group  symmetry  operations,  but  is  invariant  only  under  the  operations  of  the  ion-site  symmetry 
Gj,.  Therefore,  the  phonon  wave  functions  which  transform  according  to  the  irreducible  repre- 
sentations of  G are  combined  into  linear  combinations  which  form  a basis  for  the  irreducible 
representation  of  G^.  In  practice,  this  means  that  the  representations  arrived  at  by  using  the 
full  space-group  G are  no  longer  irreducible  in  the  less  symmetric  point  group  G^  and  may  be 
further  reduced,  analogous  to  splitting  of  residual  degeneracy  due  to  a lower  symmetry  distor- 
tion. If  this  is  applied  to  the  lattice  normal  coordinates  , the  following  expression  for 
Vgp  is  obtained : 

% = S v'  r . (16) 

kTr 

The  sum  extends  over  all  distinct  stars  of_k  in  the  BZ,  over  all  irreducible  representations 
r of  G^  which  appear  in  the  reduction  G ->■  G^,  and  over  all  rows  of  the  representation.  The 

electronic  operator  Vj^p'"  now  belongs  to  the  same  T as  Q|^p  for  a real  representation  and  to  T* 

f 

for  a complex  one.  Knowing  the  transformation  properties  of  the  Vj^p'^’s  will  lead  to  the  selec- 
tion rules  for  phonon  transitions  (whether  a phonon  can  interact  with  the  RE  ion). 

For  YAG,  the  reduction  of  the  Oj^  representations  in  the  site  symmetry  group  I)-, 
results  in  the  following:'^ 

A]  A 
A 1 B ] 

E:  .A  I + B ] 

T2  -*  A,  + B:,  -I-  B. 

T|  - B:(  -I-  B^  + Bp"  I r I 

riuis,  the  branches  shown  in  Equation  (9)  become  240  non-degenerate  vibrations 
:s  Aj,  + 3!  A^,  + 27B,g+33  Bp,  + 28  B.g  + 34  B.^, 

+ 28B3g  + 34B3„  ilM 

(tl.e  designation  of  parity  above  has  no  real  meaning  in  a crystal  without  inversion  symmetn, 
and  is  left  mereh  us  a reminder  of  approximate  parity). 

NONR ADIATIVE  TRANSITION  RATE 

With  the  above  specification  of  normal  modes  of  the  crystal,  one  can  look  at  the 
(|uanlum  mechanics  of  a direct  phonon  transition.  The  transition  probability  per  unit  time 
IS  (implicitly  averaging  over  direction  and  polari/ation  of  the  phonons) 


w 


-IT 


Wj^f  = =r  D(OI<nv,.  ii>h 


(i‘)) 


where  I)(  I ) is  the  density  of  states  defined  such  that  D(f)  dEf  is  the  number  of  such  states  in 
the  energy  range  dE,-.  Under  the  assumption  of  a low  mixing  of  vibrational  and  electronic  wave 
functions,  one  can  write  the  initial  state 


i>  - I >;'i>  IlkUr  l^kE  "kl 


, > 


CO) 


and  separate  the  matrix  element  into  two  parts.  If  we  assume  a phonon  is  absorbed  and  the 
electronic  states  go  from  I C' j > “*■  I {■  >.  then  we  can  write 


f > = I C'f>  Ilkl'r  I ^kl'  ”kr>  ^ I ^kT'  ‘"kT'  ■ ' 


Cl ) 


where  the  absorbed  phonon  in  mode  energy  ot  h ccp  - I f > - I i >• 

The  matrix  element  < f I V^p  | i > becomes 

<f|V  |i>=  V 


f'kT'r'IQkT'  n|,-p->  X lOk'T''  - I ) > 


K-) 


'^kTV  I ^k'l”  ”kT'  > 
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As  indicated  before,  Qj^p  can  be  written  in  terms  of  boson  annihilation-creation  operators  as 

. / \l/2 


=\- 


-^kr/ 


'^kl'r  ‘‘-kPr*  • 


C3) 


so  that  cciuation  (22)  can  be  reduced  to  (from  orthogonali/ation ) 


<fiV,p|i>  = <v,-|v^p'  IVi>(- 


kP 


h n 


kP 


1/2 


(24) 
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(omitting  the  primes  on  k,  P.  and  r)  1 he  iransition  probabiliiv  per  unit  time  becomes. 
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I he  motivation  for  obtaining  exact  normal  modes  and  for  deternuning  preciseU  the 
traiislormation  characteristics  lor  (he  ' cle.ir.  It  the  direct  prodiiLt  ot  the  ii reiliicilMc 

repicsenlalion  for  C/|  and  u,  iloes  not  contain  the  irreducible  representation  of  ( hence, 

the  irreducible  representation  ot  (,)|,.p  in  the  point  group  (i,.)  then  the  transition  is  not  allowed 
riuis.  by  knowing  the  transformation  properties  ot  the  initial  and  linal  electronic  wave  lunc- 
tions,  and  the  transform. ition  pri'perties  ot  the  normal  movies  ot  vibr.ition.  <vne  c.in.  via  group 
theory,  vietermme  whether  a particular  phomm  c.m  interact  with  the  Rl  ion.  eg.  whether  it  can 
cause  broadening. 
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With  the  case  of  Nd'*'^  in  low  symmetry  VAG,  the  formal  approach  proves  discon- 
certing. In  the  D-.  point  group  symmetry  of  the  Y site,  it  was  determined  that  all  electronic 
wave  functions  transform  as  E-type  irreducible  representations.  Thus,  one  must  reduce 
E X E in  D-),  which  yields,  recalling  figure  4, 


E X E = A + B3  + B + B I . 


(26) 


So,  for  the  low  site  symmetry  of  D-,,  it  is  clear  from  equation  ( 26)  that  ^ 240  nondegenerate 
modes  may  couple  to  the  Nd'''^  ioifand,  in  theory,  all  240  may  stimulate  nonradiative  decays, 
('■roup  theory  does  not  specify  the  magnitude  of  the  interaction,  however.  Obviously,  the 
approximate  odd  parity  modes  will,  in  general,  have  smaller  values  for  < f | vj^p'"  | i > than 

modes  with  approximately  even  parity.  Also,  the  nearest  neighbors  of  the  Nd'''^  ion  will  inter- 
act the  most,  ie,  the  v'  are  different  for  various  normal  modes.  With  the  formalism  well 

established,  one  may  begin  to  analyze  the  specific  example  of  Nd  in  Y AG  using  various  degrees 
of  approximation. 


INTERNAL  PHONONS 


first,  it  is  expedient  to  assume  that  the  lowest  order  interactions  are  those  with  the 
eight  nearest  neighbor  oxygen  ions.  The  .XYg  type  molecule  is  assumed  isolated  from  the 
lattice  in  that  the  molecule  has  its  own  set  of  normal  modes  of  vibration.  With  ionic  radii  for 
0~-,  Nd'*'^,  Y"''^,  and  Al'''-^  being  1 .40  .A,  1 .23  A,  0.93  .A  and  0.50  A respectively. this  is  ^ 
easy  to  accept.  What  is  planned  is  to  consider  the  eight  0~~  ions  coordinated  around  the  Nd  ' 
ion  to  be  in  cubic  0|.,  symmetry,  find  the  modes  of  vibration  of  this  isolated  complex  and  then 
reduce  the  internal  modes  of  vibration  in  Di  symmetry ; retaining  the  approximate  state  of 
parity. 

These  internal  normal  modes  of  vibration  are  not  lattice  waves,  though.  Tlie  Q s 
arrived  at  in  equation  (25>  will  stimulate  the  localized  modes  to  oscillate  and  the  localized 
modes  will  affect  the  transition  (in  our  approximation).  Thus,  if  Oj,;  is  the  normal  mode 
of  the  molecular  complex,  one  can  write  the  electron-phonon  interaction  as 


2 'k 

K 


(2") 


where  the  approximation  symbol  is  intended  as  a reminder  that  only  nearest  neighbors  are 
interacting  and  that  we  have  assumed  the  type  molecule  to  be  isolated.  Since  the  lattice 
vibrations  are  complete,  the  internal  normal  mode  can  be  expanded  linearly  in  tenns 
of  lattice  vibrations; 


‘‘K.kfr'^kr 
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where  the  aj^  ’s  are  the  mixing  factors.  This  results  in  a minor  complication  of  equation 
(24),  in  that  the  lattice  nomial  modes  must  mix  to  give  an  internal  mode;  one  no  longer  has 
the  simple  orthogonalization.  The  result  is 
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\Vj^,  = 7rD(w,.p) 
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The  above  expression,  for  the  transition  rate  of  the  direct-phonon  absorption,  can  be 
interpreted  as  follows:  A lattice  phonon  in  mode  kfr,  and  energy  h coj^p,  interacts  with  the 
molecular  complex  XYg.  This  causes  displacement  of  the  nearest  neighbors  to  the  Rli  ion  which 
have  been  described  in  terms  of  internal  normal  modes  of  the  complex  (the  amount  ot  inter- 
action being  given  in  terms  of  the  a],^  |,.p|.  coetticients).  This  variation  in  the  crystalline  tield 
induces  the  transition  of  the  electronic  state  and  the  an’.ihilation  of  a phonon  in  mode  kPr. 

The  mixing  coefficients  aj,^  |^p^  will  thus  become  very  mode-  (ie.  frequency)  dependent. 


Not  only  will  l)(cc|^p)  and  n(co|,.p)  measure  how  many  phonons  are  present  to  induce  transi- 


tions, but  the  a's  will  measure  how  much  those  particular  lattice  vibrations  attect  the  nearest 
neighbors  of  the  RE:  ion,  while  the  vj^'s  determine  how  much  the  nearest  neighbors  atleci  the 
RH  ion  (induce  transitions).  Therefore,  the  a coefficients  conceal  peaks  in  the  W(co)  curve;  a 
large  peak  is  expected  whenever  a lattice  phonon  highly  couples  into  the  internal  mode  (and 
this  is  when  the j_coefficient,  and  thus  the  (Jj^  is  large).  This  is  another  reason  tor  the  concern 
over  internal  modes.  Besides  requiring  the  transformation  properties  of  to  find  phonon 
selection  rules,  the  number  of  possible  internal  modes  affecting  the  RE-  ion  will  be  an  indica- 
tion of  the  complexity  of  the  W(oj)  curve. 

I ven  more  concealed  is  the  added  frequency  dependence  in  \V(cc)  due  to  the  poleiiiial 
gradient  term  vj^.  The  origin  of  this  dependence  is  seen  if  one  recalls  that  the  Stark  splitting 
of  each  J manifold  arises  from  matrix  elements  of  the  crystalline  potential  energy  H^..  Since 
the  coefficients  of  IE|_.  determine  the  Stark  splitting,  and  the  amount  ot  splitting  determines 
the  phonon  frequency  to.  there  is  a relation  between  (and  thus  the  derivatives  vj,^)  and  to. 

A generalization  about  the  frequency  dependence  of  < ' vj^  | i,l/j  > is  very  ditlicult. 

I lie  first  step  in  this  nearest-neighbor  solution  is  to  find  the  modes  ot  vibration  ot  this 
isolated  N't' v;  molecular-like  system.''  The  nine  ions  have  a total  of  2^  degrees  of  treedom 
and,  bv  reducing  the  reducible  representation  of  ionic  displacements  in  the  ()|j  point  group, 
we  find  the  follow ing  symmetries: ' 


Ai._.  + 1 + fi^  + 2Ts  + 


+ 1 +3T,„  + T,„. 


Jir 


(2‘)| 


If  one  I'l^i  and  one  r|„  representation  is  removed  to  account  for  transl.uiiinal  and  rotational 
ileerees  of  freedom,  the  internal,  localized  modes  remain: 


A I „ + I - I s„  + - A s^i  -I-  1:^1  + - I 1 u + T 2ir 


I .'()  I 


In  a similar  manner  by  which  eiiualion  ( I .s ) was  derived,  the  above  representations  can  be 
reduced  in  the  actual  point  group  Ds.  The  irreducible  representations  ot  0|,  iviluced  in  the 
m.inner  shown  in  equation  I”  so  tli.it  the  following  miern.il  modes  are  possible: 


4 A,^  -t-  B|g  -t-  2Bs^  -t  21)3,^  -(■  3A,|  + B | (,  + -'B2u  + 3B3^|. 
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The  distinction  with  respect  to  parity  is  used  only  as  an  approximation  (formally  we  can  say 
that  we  are  reducing  m the  group  and  consider  the  D->  site  symmetry  as  a slight  distortion 
to  D '>[,)•  Thus,  there  are  2 1 possible  vibrations  of  the  complex;  9 even  and  1 2 odd.  So 
in  the  matrix  elements  of  < 1 vj^^  | v!/j  >,  one  may  expect  9 large  values  corresponding  to  the 

9 approximately  even  parity  nondegenerate  modes. 


RAM.VN-LlkE  INTERACTIONS 

The  formal  derivation  of  the  direct-phonon  transition  rate  illustrates  tlie  complexities 
which  result  in  a crystal  of  low-site  symmetry.  This  detailed  derivation  will  not  be  followed 
in  the  even  more  complicated  derivation  of  the  Raman  broadening  term.  This  is  done  for  two 
reasons.  First,  the  inherent  formalism  of  a second-order  process  introduces  complications 
which  obscure  the  desired  result  by  requiring  the  determination  of  parameters  impossible  to 
measure  and  extremely  difficult  to  calculate.  Second,  a formal  derivation  of  the  Raman 
broadening  mechanism  is  not  necessary  to  determine  the  parameters  in  which  we  are  inter- 
ested. In  the  final  conclusion  we  will  have  to  rely  upon  experiment  anyway  to  normalize  the 
result  and  this  implicitly  calculates  all  matrix  elements  and  makes  all  the  necessary  summations. 

Just  as  in  the  direct  process,  we  desire  to  find  a transition  rate;  this  time  for  the  two- 
phonon  process  exchanging  states  between  two  energy  levels  of  a single  RE  ion  situated  in  a 
perfect  crystal.  Schematically,  we  have  the  situation  shown  in  figure  9.  For  an  upward  transi- 
tion, a phonon  of  frequency  w | is  absorbed  and  a phonon  of  frequency  cui  is  emitted.  Since 
the  transition  occurs  very  fast  (via  the  virtual  intermediate  state).  Heisenberg  broadening 
results.  As  before,  the  transition  is  described  by  an  interaction  Hamiltonian;  w'hich  means  we 
must  expand  the  crystalline  potential  as  seen  by  the  RE  ion  in  displacements  of  the  surrounding 
ions,  this  time  to  second  order. 


•VIFtTUALOH  INTERMEDIATE  STATE 


figure  'I,  Two  schematic  represenlalions  ol  the  upward  and  downward  Raman  non- 
radiaiive  energe  lecel  exchange  mechanisms 
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V = Vn  + V*  * ^ + + (32) 

ep  0 ep  ep  • 

Vq  is  still  the  static  Stark  term  and  is  the  first-order  term  used  to  calculate  (direct). 


v(  I)  = ^ 

ep  ^ 


''ep  ^ Xj  (n;o;)  Aj  (n;o!)  V I Q, 

n a i = 1,2,3 


where  Xj  (n;a)  is  the  i^*^  direction  of  the  displacement  of  the  ion  in  the  n"'  unit  cell, 

A j (n;a)V  is  the  i^^  component  of  the  derivative  of  V taken  at  the  ion  of  the  n**^  unit  cel 
and  it  is  evaluated  at  the  Ri:  ion  (since  this  is  where  the  effect  happens).  The  second-order 

expression  tor  is 

v[.“^  = -V  ^ ^ I Xj  (n;a)  Xj  (n_;a)  1 Vj  (n;o;)  Vj  (n_;£)l  VMq,  (34) 

nai  nod  [ ~ ~ ’ 

where  the  last  expression  is  the  two-component  derivative  of  V,  namely 


dxj  (n;a)  9xj  (ti^a) 


Therefore,  in  the  same  way  equation  ( 14)  was  derived,  we  can  write  Vg"*  as 


^ep  ^ ^ ''k7;  k^  ^ky  ^k7  • 

k7r  k7r  — 


- rr  y 

'k7;  k7  ''  iL-i  Aw 


(Nm^j.)  (k;  a) 


Cj,^^  (k:  a)  exp  (i  k • ( R,,|  + ))  Vj  (n; a)  Vj  ( u:  «)V  I q j . (3o) 

•Again,  R,,,  is  the  position  of  the  n"’  unit  cell  and  the  Cj.^^  (k:a)  coefficients  measure  how 
much  a normal  mode  affects  the  particular  displacement  Xj  (n;a). 

Therefore,  we  may  pet  Raman-like  transitions  two  ways:  lirst  order  in  or  second 

oriler  v[,' * . It  is  here  that  the  formalism  becomes  cumbersome  and  it  is  here  that  tljc  dcioa- 
tion  will  be  streamlined.  F rom  equation  ( 18)  we  have  (neplectinp  multiplicative  consi.mi^ 

W'K  ot  |<  f I v[,“  * I i > I-  D(i’). 


Here  D(f ) is  the  total  density  of  final  states.  F'or  an  upward  Raman  transition  we  expeet 
initial  and  final  states  which  look  like 


Ii>"  I '^i>!lkf'rl^kr  "kr> 

I I > = I > ll|,p^  I Q|.[,  n,.,,  > I - 1 )>  I (n|^..p„  + 1 I > 

for  the  process  shown  in  figure  When  the  Q|^p  's  are  written  in  terms  of  boson  annihila- 
tion-creation operators,  we  get  terms  like 

+ 

ak"p”r"  • 


which  "connect"  | i > and  | f >.  I herefore,  Raman  transition  rates  are  derived  which  look  like 


n (cu|^i'l  (n  (W|^'p'l  + 1 I 


W oc 
U 


ojkp  wp'r' 


^ 'kl':  kT' ^ *^f  ^ I 0(0)^j' • t.  (3X| 


rearranging  primes. 

As  stated  earlier,  the  purpose  of  this  calculation  is  to  find  the  temperature  dependence 
of  the  linewidth.  In  order  to  accomplish  this  rigorously,  one  must  go  through  the  same  internal- 
model  analysis  as  that  done  for  W*^.  As  will  be  shown  later,  the  nature  of  the  crystal  reduces 
the  practicality  of  carrying  out  such  calculations.  Instead,  a simple  approximation  is  to  be  used 
to  obtain  a good  fit  to  the  experimental  curve,  with  enough  versatility  retained  in  the  formalism 
to  allow  prediction  ot  other  ion-host  combinations.  As  was  emphasized  in  the  direct  case,  the 
detailed  stnicture  of  the  l)(co)  curve  tends  to  smooth  out  when  sums  are  taken  and  the  real 
density  of  states  is  used  and  especially  if  lattice  anhannonicity  is  considered.  In  the  Raman 
case,  the  integration  over  all  phonon  modes  inducing  transitions  will  introduce  even  more 
smoothing.  For  this  reason  we  can  consider  the  Raman  transition  rate  as  simply 

w|^  = (■  n(co  I I (ntcc-i  I + 1 1 fXcj  j ) DIcjt  > 


('  1 n(  CO  ] + 1)1  n(  CO  -1 ) D(  CO  ] ) l)(  (o  -I ). 


( ) 


where  in  eciuiiibrium  W,,  = W,,  ^.-AF./kT  p-  .,jng)e  multiplicative  parameter  taking  into 
account  all  matrix  elements,  summations,  internal  mode  mixing,  and  the  like.  In  practice, 
this  will  he  determined  experimentally  by  normalizing  the  linewidth  to  experimental  results 
at  one  temperature, 

Ihe  use  of  in  second  order  results  in  a similar  expression  (with  different  matrix 
elements).  Here  we  use 


^ , 

m kl'r  k'rV 


< f 1 V[.p*  I m ><  m I " I i > 1- 
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F,) 


I)(Fp-F,) 


(40) 


where  m is  a member  of  the  complete  set  of  wave  functions  comprising  the  intermediate 
states.  Since  the  difference  between  this  expression  and  the  one  used  in  first  order  only 
changes  ('  in  equation  (39).  we  will  still  use  (39)  as  the  expression  for  the  Raman  transition 
rate. 


INTERACTING  PHONONS 

INTERNAL  MOOES 

In  order  to  find  the  thermal  dependence  of  linewidth,  the  spectrum  of  those  phonons 
affecting  the  active  ion  must  be  determined.  As  was  shown,  the  existence  of  a lattice  phonon 
is  not  a guarantee  it  will  couple  into  and  attect  the  active  ion,  inducing  a non-rad iative  transi- 
tion. Acoustic  phonons  have  high  densities  ot  states  and  large  occupation  numbers,  yet  they 
correspond  to  movement  of  the  entire  unit  ceil  and.  therefore,  small  relative  motion  to  the 
active  ion  to  induce  transitions.  The  determination  of  a tractable  model  for  interacting  phonons 
is  what  is  sought  in  this  section.  If  the  range  of  strongly  interacting  internal  modes  is  deter- 
mined, the  predominant  broadening  mechanism  can  be  found  (whether  direct  or  Raman),  the 
transition  rate  can  be  summed  for  all  effective  modes,  and  a simple,  compact  linewidth  formalism 
can  be  found.  The  degree  of  experimental  fit  of  the  final  linewidth  theory’  will  justify  the 
approximations  made. 

The  Debye  theory  for  acoustic  phonons,  when  applied  to  the  system  of  a rare-c-arlh 
impurity  ion  located  in  a simple  crystal,  results  in  a linewidth  theory  well  illustrated  in  articles 
such  as  Yen.  Scott  and  Schawlow.-  This  approximation  allows  all  lattice  phonons  to  couple 
into  the  active  ion,  gives  no  clues  as  to  which  broadening  mechanisms  predominate  and  results 
in  non-analytical  solutions.  In  addition,  use  of  the  Debye  theory  requires  many  curse  fitting 
parameters  to  fit  experimental  data,  both  a multiplicative  parameter  as  well  as  a special  Debye 
cut-oil  lrec|uency  (unrelated  to  simple  bulk  thermodynamic  quantities)  for  each  set  of  levels. 

I his  theory  then  does  not  predict  the  functional  dependence  of  linewidth  with  temperature,  it 
merely  can  be  made  to  fit  the  experiment.  This  is  often  at  the  expense  of  absurd  iniplic.itioiis 
with  respect  to  physical  parameters  such  as  hulk  modulus,  thermal  conductivity . heat  capacity  . 
aiul  the  like.'^  As  will  be  shown  later,  only  the  insensitivity  of  the  thermal  dependence  of  line- 
width  with  the  structure  of  the  density  of  states  accounts  for  the  success  in  litting  the  ilieoi- 
etica!  model  to  the  experimental  results. 

l et  us  look  at  another  approach  to  the  linewidth  problem,  employing  localized  moiles. 
Consider  again  the  results  of  equation  ( 25 ) employing  internal  modes  in  the  single  phonon  case. 
Here  all  have  three  parameters.  D(cO]..p).  aj.^  and  < !.:/(•  | vj.^'  ' Wj  > - all  of  which  measure 
the  effect  ol  a lattice  phonon  ol  lrec]uency  cc|..p  on  the  active  ion.  Recall  that  Dicoppi  is  the 
ilensity  of  sttites  of  the  hit t ice  phonon  of  frequency  utpp.  a|p  pp^.  is  a me.isure  of  how  tiuicli 
the  lattice  pluinon  altects  the  localized  mode  K.  and  vjp'  determines  how  much  the  internal 
mode  atlects  the  RE  ion  and  induces  nonradiative  transitions).  Here,  for  simplicity . all  three 
parameters  are  combined  tiiul  ctilled  the  effective  density  of  states.  We  have  merely 

W'li^  = D'  (CO)  (iKw)  + I ) 

W = D'  ( CO)  II  (to), 

II 


with  similar  results  for  the  Raman  case  (equation  (39)).  D'  (to)  is  the  spectrum  ot  the  ettective 
internal  phonons  and  is  an  average  over  ^lattice  phonons  reOecting  how  the  lattice  phonons 
at  (requency  oj  affect  the  rare-earth  ion.  Now  the  non-radiative  transition  rate  is  sensitive  to 
only  these  active  internal  modes  rather  than  all  lattice  phonons. 

What  does  D'  (w)  look  like?  Certainly  it  will  not  be  a continuum  like  the  Debye 
theory.  In  first  order,  it  will  most  closely  resemble  the  Hinstein  theory,  discrete  normal  Ireq- 
uencies  for  each  internal  normal  mode  coupling  into  the  active  ion.  Figure  10  illustrates  the 
difference.  Fundamentally,  the  continuum  described  by  the  Debye  theory  is  replaced  by  a 
spectrum  reflecting  the  active  internal  normal  modes;  exhibiting  a minimum  as  well  as  a maxi- 
mum frequency  of  vibration  of  the  RE  ion,  nearest-neighbor  complex.  This  wdl  have  imme- 
diate impact  on  the  direct  phonon  transition  rate,  as  low  frequency  phonons  resonant  with 
energy  level  separations  may  not  exist  in  D'  (cc).  One  can  also  predict  that  there  is  width  to 
the  lines  in  D’  (cu)  due  primarily  to  anharmonic  effects  which  will  cause  overlap  of  peaks. 
Therefore,  all  that  is  to  be  determined  to  find  linewidth  thermal  variations  are  the  frequencies 
of  the  normal  modes  and  the  relative  magnitude.  As  will  be  seen  later,  more  important  than 
the  precise  values  of  the  vibration  frequencies  will  he  the  range  of  active  interacting  phonons  in 
D'  (cu).  With  D'  (cu)  modeled,  the  thermal  variations  of  the  linewidth  can  be  solved,  in  this 

case,  analytically. 
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Figure  10.  Four  models  for  the  density  ot  luiiiec  and  internal  phonons. 
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VIBRONIC  TRANSITIONS 

One  method  of  obtaining  insights  as  to  the  structure  of  the  effective  phonon  spectra 
is  to  look  at  the  vibrational  sidebands  accompanying  the  electric  dipole  transitions  in  RF-doped 
crystals.  Since  direct  vibrational  data  are  difficult  or  impossible  to  obtain,  this  is  extremely 
valuable.  Here  the  vibronic  data  will  reflect  the  range  of  interacting  modes  affecting  the  impurity 
RK  ion  and  will  establish  the  degree  to  which  these  modes  are  indeed  "intemal.”  By  examining 
the  structure  and  number  of  phonon-assisted  sidebands,  and  by  recognizing  that  vibronics 
are  also  a nearest-neighbor  phenomena,  one  should  be  able  to  establi.sh,  at  least,  the  frequencies 
of  the  internal  highly  coupled  modes  in  the  crystal.  Since  these  data  are  usually  taken  at  low 
temperatures  and  with  high  dopant  impurity  concentrations  to  facilitate  the  experiment,  the  ap- 
plication of  these  results  to  the  nonradiative  energy-level  exchange  problem  will  not  give  the 
exact  shape  of  D’  (co). 

Let  us  find  a simple  expression  for  the  intensity  of  the  vibronic  emission  or  absorption. 

As  before,  for  the  static  lattice  transition  in  a crystal  with  inversion  symmetry,  electric-dipole  j 

transitions  are  allowed  only  if  there  is  an  odd  parity  component  of  the  interaction  Hamiltonian  j 

to  admix  opposite  parity  states.  The  odd  parity  lattice  vibrations  will  now  cause  the  admixing. 

Consider  the  electric-dipole  transition  matrix  element 

<^,l-cH^g> 

between  the  upper  state  | > and  the  lower  state  1 vl/g  >.  If  no  perturbation  of  these  states 

exist,  the  matrix  element  will  be  zero  since  both  levels  are  derived  from  the  same  4f'^  config- 
uration and  have  the  same  parity.  The  application  of  the  electron-phonon  interaction  Hamil-  ! 

Ionian  to  the  wave  functions  will  cause  admixing  with  excited  opposite-parity  wave 
functions  from  the  4f'^“'  5d  configuration,  say,  and  permit  transitions.  Consider  the  unper- 
turbed upper  state 

I i;'u>Ilki'rlOk['»kr>-  || 

Since  we  know  that  V^.p  in  first  order  cannot  connect  states  separated  by  more  than  one 
quantum  number,  we  can  write  the  perturbed  upper  state  containing  the  admixed  state 
I > due  to  one  phonon 

lx„>=  I > >>krr  l^^kT  ”kr>  + 

- ||  I II' I n,^,.  > 1 n,^.p-  - 1>|  * V,p  | 

NC'i,'>llkr,.iC»U-  n,,p>l 


"kr 


/ 

^kr  ”kr>  ' ^\’r’  "kT'  - ' >1 


Here,  and  are  purely  electronic  wave  functions  with  the  sum  /3  extending  over 
all  excited,  opposite-parity  electronic  states.  Again,  the  lower  state  becomes 


IXg>  = \'Pg>  ' ^kT'  ^”kT'  ” 

-1^{[  I I'kPr  10|,.p  * '^ep 

* ^'^g^'^'kl>^^kr"kr^'^kT'*”kT'~ 

E(3  - Eg 

Here  a simplification  is  made.  Since  most  vibronic  data  are  taken  at  low  tempera- 
tures (to  avoid  second  order  effects  such  as  anharmonics)  the  limit  of  T = 0°K  is  assumed. 
Thus,  the  occupation  numbers  for  ail  modes  tend  to  zero.  Therefore,  using  the  first  order 
expression  for  given  in  equation  (14),  orthogonalization  of  the  wave  functions,  and 
retaining  only  first  order  terms,  we  obtain 

<Xu  I - e r I Xg  > = 


l'^|3>^krrlQkr'’kr>  ■ 


Ejj  - Eg 


/ 

< I V|,p  I X *^i3  I ■ ^ ^ ' *^g^ 


(43) 


I his  corresponds  to  a transition  of  the  RE  ion  from  state  g to  state  u plus  the  creation  of  one 
phonon  m mode  (ki'r).  Here  we  see  that  1 1/^^>  must  indeed  be  derived  from  an  opposite 
purity  configuration  to  remain  non-zero. 

The  transition  probability  per  unit  time  is  given  by  the  familiar  relation 


c r IXp>  P<E)  . 


where  p(  I ) is  the  density  of  final  states  composed  of  the  product  of  the  lattice  phonon  den- 
sity of  states  l)(tu)  and  the  electronic  states  g(E  - h).  Since  the  energy  spread  of  the  elec- 
tronic states  is  much  smaller  than  the  range  of  the  phonon  states,  g(E  - hcu)  = 6(E  - hcu  - 
cU|^p)  and  the  vibronic  spectrum  retlects  the  structure  of  D(oj)  (figure  1 1).  The  precise 
structure  o(  IXoj)  is  diUiciilt  to  determine,  though,  becau.se  of  the  electronic  matrix  elements 
in  equation  (43)  and  because  of  the  difference  in  the  matrix  elements  involving  'kT 
example,  since  the  is  opposite  parity  to  and  v!/g,  only  the  odd  parity  vibrations  can 
have  non-zero  matrix  elements.  Recall  that  the  expression  for  the  single  phonon  nonradiative 
transition  depends  on  the  even  parity  normal  modes.  The  group  theoretical  selection  rules 
are  also  similar  to  those  for  phonon  nonradiative  transitions.  This  states  that 


w 


30 


a(CM" 


NOTE:  A,  B,  C,  D,  E AND  F ARE  WEAK  ELECTRONIC  TRANSITIONS 
AND  ARE  TO  BE  REMOVED  FROM  THE  VIBRONIC  SPECTRA 


Figure  1 1 . Vibronic  spectra  at  low  temperatures  for  Yb^  A1-,  (AI04)3 
indicating  effective  phonons. 
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tioes  not  contiiin  an  irrotlticihle  representation  of  the  tratisitioti  ntotnent  r when  reduced  iti 
the  site  symmetry  (where  1'^,  atid  1'^,  are  the  irreditcihle  representations  ol  the  upper  and 
ground  stares  respectively,  and  is  the  irreducible  representation  of  the  normal  mode 

of  Cij,).  In  crystals  such  as  Y.Afi.  which  tio  not  contain  inversion  symmetry,  parity  is 

not  a good  tiuantum  number  and  nornuil  moilcs  are  cornposerl  of  both  even  and  mid  p.irity 
parts.  Also,  the  low  D-.  site  symmetry  docs  not  exclude  any  transitions  due  to  the  selection 
rules  and  all  modes  can  catise  sidebands.  Wh.it  is  significanl  is  th.it  now  the  range  of  inter- 
acting phonons  can  be  determinerl  by  extimining  the  vibronic  spectra.  Hy  counting  the 
number  of  sidebands  and  comparing  this  to  the  number  predicted  from  the  internal-mode 
tipproximation,  the  degree  of  isolation  of  the  internal  imtdes  can  also  be  imesiigated 


PRI  .SlJMtl)  INTER ACTINC;  EFFECTIVE  MODES 

In  order  to  find  the  effective  phonons  of  the  system  til  Nd'''-'  in  ^ .Ati.  the  Mbrmtic 
spectra  at  low  temperatures  can  be  examined  The  data  available'"  are  for  A'b^  Als  ( \lO4l; 
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Since  the  ditTerence  in  atomic  weight  between  and  Nd"*^^  is  1 73  vs  144,  this  data  is 
assumed  to  be  adequate,  though  not  precise.  Also,  Buchanan  states  that  using  Yb  Al  G 
instead  of  Y'''^  doped  YAG  results  in  larger  absorption  coefficients  without  seriously  per- 
turbing the  symmetries  and  so  facilitates  the  experiment.  The  results  are  shown  in  figua*  1 1 . 
Mere  we  find  approximately  20  large  sidebands  occurring  between  250  cm“*  and  570  cm”' 
and  a small  group  of  approximately  nine  smaller  ones  from  75  cm”'  to  225  cm”'.  Recalling 
the  results  of  the  group  theoretical  reduction  of  the  XY^  complex,  we  expect  the  internal 
normal  modes  to  consist  of  nine  approximately  even-parity  and  12  approximatel>  odd-parity 
vibrations,  a total  of  21  close-coupled  modes  corresponding  to  the  XYg  molecular  system  csf 
Nd"*"-^  and  0“-.  The  conventional  normal-mode  analysis  of  the  80  ions  in  the  unit  cell  predicts 
240  possible  modes  seen  by  the  Nd'*'^  ion  in  Dt  site  symmetry.  The  data  using  Yb  Al  G fit 
well  with  the  predicted  results.  The  range  of  internal  vibrations  then,  is  confined  to  a region 
approximately  300  cm”'  wide  containing  narrow  peaks  corresponding  to  highly  coupled  modes 
(contrasting  sharply  with  the  Debye  prediction).  The  small  peaks  on  the  low  frequency  side 
are  thought  to  be  next  nearest-neighbor  interactions  of  low  energy  between  the  RE  ion  and 
the  Al'*'-^  ion.  Since  the  atomic  weight  of  Al  is  27  and  that  for  0”-  is  16,  we  should  expect  a 
lower  vibration  spectrum.  Also,  the  larger  inter-atomic  spacing  to  the  Al'*'^  ion  will  result  in 
lower  fre<)uency  mcxles  as  well  as  smaller  coupling  (as  seen  by  the  amplitude  of  the  sidebands). 
Therefore,  let  us  choose  to  model  the  complicated  effective  phonon  spectrum  by  a simple 
uniform  density  of  modes  ranging  from  approximately  250  cm”'  to  570  cm”'  (note  figure  1 2): 


1)  250  cm"'  <htL)<  570  cm  ' 

0 otherwise 


(44) 


EFFECTIVE  PHONON  ENERGY 

f igure  12.  t'niform  density  model  for  the  effective  modes  in  Nd’*’-^  Y.A(i. 
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(with  D a constant ).  Since  incorporation  of  the  precise  shape  (if  it  could  he  extracted  from 
the  data  available)  would  result  in  another  numerical  result,  let  us  try  the  above  approxima- 
tion and  examine  the  quality  of  fit  to  the  experiment.  After  the  next  section  we  will  examine 
the  prediction  of  thermal,  linewidth  variation  incorporating  the  above  approximation. 


DLBYt  LINHVVIDTH  THEORY 


What  is  presented  in  this  section  is  a brief  overview  of  the  Debye  theory  for  the 
thermal  dependence  of  linewidth,  best  presented  by  Yen,  Scott  and  Schawlow.^  This  is  done 
primarily  to  place  iti  proper  perspective  the  theory  to  be  presented.  Also,  this  turnishes  a 
theoretical  has^  to  the  effective-mode  theory  to  be  presented  later.  Assuming  a completely 
Debye  model  for  the  crystal,  Yen.  et  al,  begin  by  writing  the  Hamiltonian  of  the  impurity 
ion  in  a host  cr>stal  as 


'tot 


= H, 


"p^ 


where  IIq  describes  the  ion  in  the  presence  of  a crystalline  field.  Up  describes  the  phonon 
system 


*P  " “k'\  ‘‘k 


and  Hi  describes  the  interactions  between  the  ion  and  the  lattice  vibrations 


I',  + 

Cwj,  (a,.-a|,)  + 


D CO,..'  (a^^  - a,;,  Ha,.'  -a^.'  ) + . . 


Here 


'k 


= Hj  -*■  ir,’  ^ . . . 

and  a^  are  annihilation  .iiul  iicalion  o|ierat('rs  ol  acoustic  phomms  in  the  k’*'  mode. 
CO,.,  is  then  frequence , and  ('  and  I)  .ire  (he  linear  and  quadratic  crxstal  field  coupling  opera- 
tors respectively  (which  operate  on  electronic  states  of  the  ion  onlv  I.  M is  the  mass  of  the 
crystal,  and  v is  an  average  sound  vclocite  in  the  cresial.  Since  this  assumes  the  Debye  approx- 
imation. V is  assumeil  independent  ol  kand  k = cu  v.  for  the  Ivpical  energv  level  conligura- 
tion.  the  relaxation  transition  prob.ibilitx  lor  the  I'l'  level,  due  to  the  direct,  single  phonon 
process  (Wj*)  is  (using  ’t  en's  luXalion) 


W^'  = 


J7T 


\ 


s.  1 — 


:Mv-h|,<i 


CJ,I  /XiCp 1 1<  \|  ( I A| > |- 


UycOij)+l!  + X Wjj  p(cJij)|<Aj|ClAj>|-Po(Wij) 

i>  i 


= tt)  1 ^3jj  |p^,(Wjj)  +11+1  /3jj  Po(wij)[  . 

/ J < i J > 1 > 

Mere.  p(eOjj)  is  the  detailed  Dehye  density  of  states  at  cOjj.  and  Po<<aJjj)  is  the  phonon  oeeii- 
pation  number  for  phonons  of  energy  httOjj)  at  temperature  T. 

The  contribution  to  the  relaxation  process  due  to  Raman  scattering  of  phonons  b> 
the  impurity  ion  is  calculated  using  H’|  in  second  order  and  H'j'  in  first  order.  This  results  in 
an  expression  for  the  i*^’  level  given  by 


a-  (T'0n>^ 
1 h ' 


(e’^  - 1)' 


where  x = hcc  'kT  and 


KAj|C|A|.>|-|<A,,|C|Ai>|-  ■ 

^ jk  ^ki 


5 f,  V-  I < A]^  I C I Aj  > 1 “ 

T^hr)  y . 

2 TT-  p V L ^ 'Is 

where  Oj)  is  the  Debye  temperature  (which  turns  out  to  he  unrelated  to  the  values  consis- 
tent with  bulk  data  such  as  specific  heat).  The  integral  T)  has  been  tabulated  by 

/iman  and  is  proportional  to  ((?[yT)^  at  high  temperatures. 

Ihe  treatment  of  inhomogeneous  strain  broadening  and  homogeneous  multiphonon 
relaxation  broadening  is  essentially  identical  with  that  presented  earlier.  For  the  i’^^  level, 
the  total  contribution  to  the  linewidth  due  to  M]  is.  therefore,  given  by  Yen  to  be 


An,  = ^ l^ij  (iV^ii* 

j < i j < I 


+ 'r  1 p,,(w,j)  + — j^lTO,))'  ?(,Mli)T) 

i>i 


.^4 


J 


wlKTf  K||  Is  the  (essentully  ) temperature  mdepetulent  multiphonon  contribution.  For  a given 
optical  transition,  the  Lorent/ian  contributions  I'rom  each  level  add 

A/^ij  = Ai'i  + Ar-j. 

1 he  contributions  trom  the  (iaussian  shaped  strain  broadening  are  then  added  to  Ar'y. 

I he  lit  ol  this  theory  to  experiment  is  difficult.  Since  only  Ar-y  is  actually  measured, 
the  ImesMdthos-tetnperature  curve  must  be  measured  for  al]_optical  transitions  to  separate  out 
/j'l  trom  *Vj  and.  hence,  the  various  multiplicative  coefficients  (/3jj.  Qj,  and  Ky  for  each  level  I. 

In  .iddition  to  this,  an  elaborate  fitting  technitiue  must  be  used  to  determine  the  proper  value 
tor  In  general,  a multiparameter  fit  reipiiring  various  p's  and  a's  as  variables  is  employed 
anil  these  parameters  are  adjusted  until  a good  fit  is  obtained.  For  example,  in  the  case  of  the 
^l’,j  • < s ^ transition  of  in  La  F^,  Yen  was  able  to  reduce  the  number  of  coef- 

licients  trom  (•>  to  3.  belbre  beginning  the  parametric  fit  (which  depends  on  the  dramatically 
dilierent  temperature  dependencies  of  the  direct  and  Raman  processes).  The  fit  to  experiment 
IS  c\i  client 

I he  obiective  of  the  Debye  theory  is  not  to  predict  the  shape  of  linewidth-vs-temperature 
cunc  lot  Kl  ions  m cry  stals.  nor  is  it  to  fonnulate  an  analytic  expression  for  linewidth-vs- 
tcmpcr.iture,  but  r.ither  to  check  the  validity  of  the  proposed  broadening  mechanisms  them- 
scbes  idiiect  phonon,  Raman,  etc).  The  success  of  demonstrating  this  has  been  excellent. 

I he  tiadcotl  IS  that  the  theory  has  four  constraints  relating  to  its  utility  of  being  a means  of 
predicting  the  shape  of  the  hnewidth-vs-temperature  curve  with  only  a few  measured  points: 
the  theory  of  't  en  uses  the  Deby  e theory  which  is  not  valid  for  most  hosts  (where  optical 
phonons  interact  strongly  with  the  impurity  ion)  and  retjuires  complicated  curve  fits  to  deter- 
mine an  etieciive  different  from  the  results  of  bulk  measurements;  the  functional  expression 
tor  the  R.iman  contribution  is  an  integral  form  resulting  in  a nonanalytic  expression;  the  't  en 
theory  requires  the  measurement  ot'many  points  on  many  linewidth-vs-tetnperature  curves  to 
enable  the  complicated  paratnetric  fit  to  be  completed;  and  finally,  the  theory  just  does  not 
■predict"  the  slope  of  the  curves.  The  theory  to  be  presented  will  address  the  prediction  ol 
the  shape  of  the  linewidth-vs-temperature  cuive  for  rare-earth  impurity  ions  residing  in  hard 
cry  stals.  Since  data  are  available,  and  because  of  its  wide  application  as  a laser  medium.  Nd^"' 

> ,\(  I u ill  be  investigated  specifically . The  regime  of  300  K to  400°K  is  of  special  interest  as 
this  is  the  ty  pical  operating  temperature  of  high-power  Nd^'^  \ ,\(i  lasers  and  amplifiers,  l lie 
pr.ictical  utility  of  a simple  expression  for  the  thermal  dependence  of  linewidth  a-quiring  only 
the  measurement  of  the  300°K  linewidth  and  the  low-tcmpcraturc  intercept  for  each  optical 
transition  desired  is  obvious. 


FFFI  ( TIV  F MODF  I.INFVVID  IU  TIIFOm 
RAM  AN  I FFF(  T OVF  R AI  L MODFS 

1 he  incorporation  ol  all  the  preceding  quantum  electronic  formalism  into  ,i  theory  li' 
determine  ihe  thermal  variation  of  the  homogeneously  broadeneil  railiative  iransiiions  used  in 
the  operation  of  Rl-  lasers  requires  the  determination  of  the  nonradiative  energv  level  exch.mec 
rale  iH'lween  all  sets  ol  levels  via  all  mechanisms.  This  rei]uires  that  the  fundamental  depen- 
dence of  linewidi h-vs-icmperature  incivrporating  all  effective  modes  for  each  mech.iniNin  must 
be  loiiiul.  the  detailed  set  of  energy  levels  must  be  examined  to  find  .ill  groups  ol  inieivicimg 
levels.  ,md  those  mechanisms  whicli  can  be  neglected  must  |ve  determined.  Hy  utih/mg 


3,‘s 


tile  new  approximation  lor  interacting  modes,  tor  example,  the  predominant  broadening 
mechanism  (in  Nd'*'-^  VAC.)  will  be  shown  to  be  the  Raman  type  and  an  expression  will  be 
derived  tor  the  total  nonradiative  transition  rate  due  to  the  summation  of  all  effective  modes. 
This  will  also  set  a practical  limit  for  the  maximum  separation  of  energy  levels  connected  via 
the  Raman  process.  An  expression  will  then  be  derived  for  the  elastic  Raman  scattering 
process  which  homogeneously  broadens  an  individual  level  independent  of  all  others.  This 
will  turn  out  to  be  a very  important  mechanism  and  will  greatly  simplify  the  anal\  tic  expression 
used  for  the  thermal  variation  of  linewidth. 

It  IS  very  important,  though,  to  remember  that  the  results  to  he  presented  are  not 
applicable  to  ion-host  laser  systems  in  general.  In  fact,  even  the  same  rare-earth  ion  in  another 
crystal  may  fall  into  regimes  not  accurately  described  by  these  approximations  (neither  quali- 
tatively nor  quantitatively).  The  system  of  Nd'*'^  YAG  when  used  as  a laser  medium  assumes, 
among  other  things,  a RF  ion  (of  small  spatial  extension  causing  low  mixing  of  electronic  and 
vibrational  wave  functions)  located  in  small  quantities  (reducing  ion-ion  interactions)  in  a 
single  crystal  (low  impurity  levels  to  form  perturbation  centers)  with  a large  unit  cell  (elimina- 
ting acoustic  mode  effects)  in  a ■‘hard”  crystal  (allowing  large-frequency  phonons  to  interact 
with  the  ion).  Also,  the  extremely  low  site  symmetry  without  inversion  relaxes  most  selec- 
tion rules  and  a system  with  the  impurity  ion  having  an  atomic  weight  much  larger  than  its 
neighboring  ions  simplifies  the  vibrational  structure.  .So,  while  the  theorx  to  be  presented  here 
can  probably  be  modified  to  satisfy  the  general  case  of  an  ion-host  combination,  rederivation 
from  fundamentals  may  be  necessary. 

First,  the  functional  expressions  must  be  found  for  the  Raman  broadening  term,  sum- 
ming the  effects  for  all  active  modes.  The  upward  and  downward  Raman  nonradiative  tran- 
sition rates  between  tw'o  energy  levels  separated  by  Al.  are  given  by  the  taniihar  relations  (in- 
elastic phonon  scattering) 


= C n(  w I ) I n(uj  s ) + I 1 !)'( to  | ) D'(co  s ) 


with 


W 


-K 


r n(  to  s ) 1 n(  CO  ] ) + I I 1)'(  to  ] ) 1)'(  tos ). 


w ith  to  ] - tos  = AFi  (note  figure  ^)),  If  all  processes  w ithin  the  operation  of  this  system  as  a 
laser  occur  with  much  larger  time  constants  than  the  nonradiative  exchange  rate  (^10”'  ' sec), 
then  it  may  be  assumeil  that  the  levels  are  always  in  thermodx  namic  equilibrium  and  described 
by  Holt/.man  statistics.  Therefore.  and  Vv|^are  not  independent.  Notice  that  this  approx- 
imation breaks  down  when  laser  mixle  locking  occurs  (with  pulse  widths  10“'  ~ sec)  and  when 
the  level  separation  becomes  l.irge  and  the  nonradiative  exch.mges  are  no  longer  described  b\  ,i 
firsi-ortler  effect;  eg.  by  multiphonon  effects.  I he  eiiuilibrium  of  the  41  j | s ihermal  level 
anil  the  4lq  s ground  state  is  ’wo-phonon  process  with  a relaxation  time  of  approximateb. 


0.5  X 10”*’  seeond  and,  therefo/e.  in  the  dynamic  opir.ition  of  a (.t-switched  laser  (r  10”*''  sec) 
these  levels  are  far  from  remaining  m thermodynamic  equilibrium  But  tor  all  the  firsi-order 
broadening  mechanisms 


where  AF  is  tlie  energy  level  separation  between  the  upper  and  lower  states  at  T.  If  the  equili- 
brium occupation  numbers  are  substituted  into  the  Raman  relations  above,  the  following  is 
obtained  (noting  that  n(cc)+  1 = n(u)) 


O'  - 1 a/3  - 


, , , AF^ 

D (tc-i ) I)  VCO-l  +-T") 
- h 


\V>^  = r -^-4^  I)’(cc,)I)'(rc,  —) 
u a - 1 a/3  - 1 ' ‘ ti 


tlOJ  I kl 


/3  = e-^' 


Therefore,  all  that  remains  is  to  integrate  the  above  expressions  over  all  interacting  modes  to 
determine  the  contribution  from  all  modes.  It  a detailed,  highly  structured  expression  tor 
the  effective  density  of  states  was  used,  it  is  clear  that  the  above  expressions  would  be  sensi- 
tive to  peaks  m D'tccl  separated  by  AF.  It  is  also  clear  that,  for  many  such  peaks  in  D'(cc». 
the  results  would  tend  to  average  to  those  obtained  by  employing  a uniform  effective  density 
of  states  (the  approximation  to  be  used). 

Hefore  integrating  the  above,  let  us  find  the  maximum  energy-level  separation.  Al  . 
connected  by  the  Raman  mechanism  employing  the  unifonn-density-of-states  approximation. 
Here  we  normalize  equation  (44) 


l)'(  CO  ) - 


‘^min  ^ ^ ^ ‘^'max 


0 otherwise. 


and  cOp,.,j^  = 570  cm' 


where  in  this  c.ise  cO|,,j„  = 750  cm~'  and  cOp.|.|j^  = 570  cm  . Note  that  D (co)  is  not  inte- 
grable  over  all  co’s  to  give  three  times  the  number  of  ions  per  unit  cell,  since  these  are  oni> 
the  effective  modes.  I'he  Debye  approximation  normalizes  in  this  manner  and  therefore 
considers  all  nu'des  eqti. ills  ettci.  livc’.  .Since' il  is  always  true  that  co | — cos  = AF  . tor  VV  j it 


must  be  that  to s '•  co,,,| 
example  ol  Nd^'^  ^ \<  ■ . 


in  that 


or  to  s < co„,y^  - Al  . Therefore,  for  the 


Al  ■ .770  cm  ' 750  cm  ' < co  s < 5'^0  - Al- . <4S) 

riierctore,  al  least  to  first  order  m our  approximations,  levels  separated  by  over  .770  cm”'  are 
not  connected  via  the  R.im.iii  mechanism  and  .is  Al.  -*  770  cm”'  the  quantity  ot  eflectise  modes 
interacting  with  the  levels  tends  towards  zero.  I or  upward  Raman  transitions,  .igain  to  | - tos  = 
Al  .md  since  tos  > tO|„u,.  to  | > tO|,,|,.,  ■(•  AF  anil  to  | < to„.my.  therelore 


770  ciir 


■t  750  cm  ' < to  1 < 5’7()  cm”' . 


The  niniie  (h'  values  for  cos  and  cj  j in  equations  (48  ) and  (49)  now  specify  the  limits  ol  inte- 

“ R R 

gration  in  equations  (4(i)  and  (47)  for  Wj  and  VV^  . 

As  an  example  of  how  the  structure  of  D'(co)  is  used,  consider  the  Stark  doublet 
forming  the  upper  lasing  level  in  Nd^^  YAG  — the  4F2^s(2)  and  4p2^2*  (note  figure  5). 
Since  the  closest  set  of  levels  are  the  dFj  /-.  (higher  in  energy  by  approximately  870  cm“' ). 
these  two  levels  can  be  considered  isolated  and  the  nonradiative  relaxation  rates  between  them 
will  determine  the  width  of  each  level  (neglecting  elastic  scattering  for  the  moment).  Since, 
from  our  approximation  we  find  that  D'(88  cm“‘ ) = 0,  the  direct  single-phonon  transition 
rate  is  zero.  Realistically,  is  considered  small  and  may  be  neglected  (this  will  be  dis- 
cussed later).  Therefore,  the  Raman  expression  alone  defines  the  broadening  of  these  levels. 
Using  equations  (48)  and  (49)  the  width  of  the  4F3^2^-' 


Ai^ldF^/itZ)] 


+ ^ strain 


all 

effective 

modes 


where  Equation  (46)  is  summed  for  250  cm  '<co-i<482  cm  *.  Likewise, 


Z)a^[4F3  ,'2(1)1 
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^ strain 


all 
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modes 


where  equation  (47)  is  summed  for  3,48  cm”'  < coj  < 570  cm”' . Ai'strain  inhomo- 
geneous residual  linewidth  due  to  microscopic  strain  within  the  crystal.  Properly,  the  Lorentzian 
shaped.  Raman  broadened  line  should  not  be  added  to  the  Gaussian  strain  line  as  linewidth  is 
defined  differently  for  each.  The  error  will  become  more  pronounced  at  low  temperatures 
where  VV'^  0.  Also,  it  should  be  noted  that  the  two  residual  strain  widths  need  not  be  equal. 

The  integration  of  equation  44  and  45  is  straightforward.  Using  equations  (48)  and 
(49).  the  following  is  obtained: 
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w here  a = e'’^  '' ' . /3  = p'  = 

incorporated  into  the  limits.  Since  W 
45).  let  us  just  integrate  , 
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is  related  to  by  the  Holl/man  ei|uatii>n  lequ.ition 


If 


f,,  I'kT  . . , 

a = e‘  ^ ‘ then  da  = — dw, 

kT 


so 


u h 


/ 


1 I 


da, 


a ' 1 ajS  - \ 

/itii  primes  eliminated  and  (i  = This  expression  becomes 


u h 


/ 


■'max 


da 


^min  + /3a-  + : 1 - 


wliich  is  in  the  form 

/dx 

A + 2Bx  + Cx- 


with 


(50) 


The  following  results  depend  upon 


AC<  B-. 


/3< 


(^) 


A 

This  is  seen  to  always  hold  as  the  above  reduces  to  ]3-  2;3  - 1 and  since  0 < /i  < 1 (Al 

this  conilition  is  always  satisfied.  Therefore  the  integration  results  in 
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So  substituting  for  A.  If  atul  C iti  eijuatioti  I 50)  and  sitnplify  ing.  the  above  expressiiMi  bec<Mnes 
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If  the  limits  are  included  and  the  expressions  fora  and  (3  substituted,  the  following  final 
result  is  obtained: 
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VV  , , when  integrated  over  the  proper  limits,  is  given  by  e 
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times  W 
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(53) 


1 xpressions  (52)  and  (53)  are  analytic  expressions  requiring  only  the  determination  of  the 
multiplicative  parameter  C (through  fit  to  a single  temperature)  to  determine  the  functional 
dependence  of  the  Raman  transition  rate  with  temperature. 


RAMAN  ELASTIC  PHONON  SCATTERINCi  APPROXIMATION 

The  elastic  scattering  of  phonons  from  the  impurity  RI-  ion  also  indiiees  a homogeneous 
broadening  of  that  level.  Rather  than  being  an  independent  broadening  mechanism,  this  process 
may  be  thought  of  as  lifetime  broadening  due  to  the  nonradiative  transition  rate  between  two 
levels  whose  energy  separation  is  zero.  In  fact,  since  all  levels  in  Nd^^  A'ACi  are  twofold  Kramer 
ilegenerate.  elastic  scattering  ma\  be  thought  of  as  inelastic  scattering  between  levels  whose 
separation  approaches  zero.  The  expressions  derived  in  equations  (54)  and  (55)  ma\  be  used  to 
determine  the  itulividual  w idth  of  a single  level  by  looking  at  the  limit  as  AE  -*  0.  Contribu- 
tions to  the  Raman  component  ot  linewidth  broadening  will  come  from  AI'  = 0 as  well  as  frcnn 
AI  ecjual  to  energy  level  separations  to  adjacent  .Stark  levels.  I inploymg  equation  (51  ).  it  is 
seen  that 


Thcrelore.  1'  Hospital’s  rule  is  employed  and 


and  that,  at  high  temperatures,  « T-,  which  agrees  with  the  limit  obtained  using  the  Debye 
approximation. 

Rather  than  requiring  both  t'orins  of  the  Raman  transition  rate,  an  approximation  will 
be  employed  which  will  greatly  simplify  the  actual  determination  of  linewidth-vs-temperature 
for  a general  set  of  Stark  levels.  This  approximation  is  to  employ  only  lor  the  Raman 
broadening  component  of  a level.  Consider  figure  I.  the  energy  level  diagram  for  a general  Rl 
ion  in  a cry  stal.  If  the  Raman  contribution  to  the  broadening  of  the  i**’  level  is  desired,  say . 
the  following  is  found: 


= ^ l7i  W*^(i)  + '^'1^ 


h> ' Ip 


where  7,  is  a multiplie.ilive  factor  related  to  the  magnitude  of  the  cimtribution  from  el.istk 

seatteriim  ( related  to  matrix  elements).  U is  the  upward  Raman  transition  rate  from  the  1''' 

" ^ K 

level  to  the  a"'  level  with  a magnitude  7^^^  (this  is  relateil  to  C in  I ij nation  ( 50)).  and  W | is  the 

downward  transition  rate  Irom  the  i’*'  level  to  lliefj’''  level  with  a magnituile  7j^,s.  I he  suins 

extemi  over  all  levels  alrove  Jj'  • ami  ail  levels  below  the  i'^^  level  with  7|^^  = 7|^’  = 0 tor 

AT  > .^dOcnr'.  Since  tlie  (luanlimiive  taclois  7^^^  and  7j^  are  impossible  to  determiiic  wiilioui 
precise  knowleilge  of  the  matrix  elements,  anil  since  their  v.ilue  is  not  expected  to  v.iry  di.isi- 
ically  with  the  wave  functions  used  within  a multiplet.  all  the  7m's  and  7,^j’s  .ire  .issumed  equ.il 

If  this  is  assumed,  a curve  of  \v|^  vs  AT  is  plotted  (tor  .1  particular  I ) and  an  estim.ite  oi  the 

^ R • 

relative  contributions  of  each  mech.imsm  can  be  obt.imed.  F'igure  1,5  shows  \s  ' I I01 

^ l\ 

'[' = .5()()  K and  T = 100  K Both  show  substantial  reductions  in  the  coniributn'n  ot  '5^^  vs 
VV*^  lor  typical  values  of  AI  say. AT  ' 100  cm“'.  .\lso.  fivr  AF.  < 1 50  cm“' . plots  of  W vs 

I and  \v|^  vs  F'  differ  by  less  than  10  percent.  Therefore,  for  Al  < about  1 50  cm~'  the  dit- 
lerence  ol  Ai'  vs  I'  using  the  two  tnethivds  is  lower  than  experimental  accuracy  ,md  lor  .M  ■ 

I 50  cm”'  the  eonlribiition  ol  mterlevel  Kviman  broadening  is  considered  small.  File  expression  lor 
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Figure  13.  The  upward  Raman  luniradiative  transition  rate  (Wj^)  vs  energy -level 
separation  AF  t'or  temperatures  of  I00°K  and  300°K. 

O 

clastic  Raman  scattering  W*'  given  in  eciuation  (54)  will  he  used  to  determine  the  temperature 
ilependetice  of  the  Raman  two-phonon  contribution  to  linewidth.  Again,  as  with  most  ap- 
proximations made,  the  tittality  of  the  predicted  fit  to  the  experitnental  results  will  indicate 
the  validity  of  the  approximations  made.  In  addition  to  this  recent  assumption,  many  approx- 
imations and  assumptions  were  made  to  derive  the  expressions  for  the  Raman  transition  rate: 
complete  local  thermodynamic  equilibrium ; a uniform  density  of  states  described  by  the  inter- 
nal mode  approximation  with  limited  ranges  (ie.  equation  (44));  and  a completely  first-order 
tiescription  of  lattice  phonons  which  neglects  anharmonics  and  other  higher-order  effects  which 
obscure  crystalline  wave  functions  and  induce  mixing  of  electronic  and  lattice  wave  tunctions. 
induce  an  appreciable  phonon  mode  broadening,  and  generally  destroy  first-order  theories  as 
ionic  displacements  become  large  due  to  elevated  temperatures. 


DIRECT  TRANSITION 

Since  the  width  of  a typical  energy  level  is  much  smaller  than  l)'(w)  (approxim.itely 
.3  cm"*  vs  .300  cm“' ).  the  expressions  for  the  direct  phonon  mechanisms  are  governed  b\ 
resonant  transitions, 

VV*^  = ( 'n(w)  l)'(w)  fidicc  - AT.) 
u 

= (''(n(cc)+l|  l)'(to)  6(hcc  - AT  1. 
d 


4: 


Tlieret'ori;, 


r 


u ,.AH/kT_ 


D'(AE/h), 


W'V  = ,,T- l)'(AE/h).  (56) 

li  ^.Ah/k1  _| 

with  l)'(Al,,  h)  given  hy  liie  norniali/.ed  equation  (44).  For  AE  < and  AE  > cOf^iax' 

''  u “ "d  ~ ^ (more  properly,  becomes  negligible).  This  greatly  reduces  the  contribution 

ol  to  level  broadening  since  the  expressions  are  governed  by  phonon  occupation  numbers 
n(w)  which  are  large  tor  small  AE/kT,  consequently,  D'(AE/h)  tends  to  render  the  upward  direct 

phonon  process  (W*^)  inetTective  for  the  regions  of  maximum  effect  - low  energy  spacings 
o _ I 

with  large  n(te)’s.  Also,  since  most  Stark  levels  are  separated  by  less  than  250  cm  , few  levels 
even  tall  in  the  range  of  l)'(to)  (and  can  couple  with  the  effective  phonon  via  the  direct  mech- 
anism). If  Ai;  > coj^in  ‘^max-  allowed,  but  n(cc)  tends  to  be  small  at  all  but 

the  highest  temperatures  and  wfj*  remains  relatively  ineffective  in  broadening  the  levels.  Notice 
also  that  for  large  temperatures,"  the  Raman  mechanism  increases  much  more  rapidly  than  the 

direct  a T-  as  compared  to  wjj*  T). 

W*?,  on  the  other  hand,  tends  toward  a constant  (of  unknown  magnitude  relative  to 
d 

other  mechan/.sms)  as  nioj)  hecomcs  .small  and  therefore  contributes  little  to  the  thermal 
variation  of  the  linewidth,  and  could  he  lumped  together  with  the  temperature  independent, 
strain  induced  linewidth.  This  temperature  independent  contribution  to  linewidth  is  then 
determined  independently.  For  the  above  reasons,  the  temperature  dependence  of  the  direct 
phonon,  nonradiative  transition  broadening  mechanism  is  considered  to  be  small  when  com- 
pared to  the  Raman  process  when  determining  the  thermal  variation  of  linewidth  in  a system 
of  Ntl'*'-^  ^’A(i.  Since  the  multiphonon  process  is  also  small,  the  thermal  broadening  of  two 
energy  levels  is  assumed  to  be  due  to  residual  strain  and  the  Raman  process  only,  fhis  becomes 
another  assumption  to  be  redeemed  by  the  ciuality  of  fit  to  the  experimental  data  (and  does 
imlectl  break  down). 

(n  passing,  one  should  note  two  injporlani  points.  First,  for  large  T's  (above  40f)'’K 
say),  this  approximation  will  begin  to  fail  as  modes  between  250  and  570  cm”'  will  become 
significantly  populated.  Therefore,  in  addition  to  the  effects  of  anharmonics  and  other  higher 
order  effects,  the  high  temperature  fit  is  expected  to  degrade  due  to  the  assumption  of  small 
lunctimial  dependence  of  W*\  Second  and  most  important,  for  some  intermediate  energy 
level  spacings,  ilic  assum|ilion  ot  a single,  uniform  tlensiiy  of  states,  with  limits  ol  250  cm”' 

.iml  ^''0  cm”'  becomes  important.  Here  it  was  assumed  that  the  next  nearest-neighbor  modes 
(approximately  one  order  of  magnitude  smaller)  from  75  cm  to  225  cm”  were  negligll^le 
I his  issumption  ni.iy  prove  invalid  it  very  large  matrix  elements  (reflected  in  C ')  occur  I'n- 
fortim.ilcK  . prciliction  when  exii.iordm.uily  large  values  ol 
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occul  iM-\liemclv  dilliciilt.  Ihercfore,  the  approximation  will  be  used  with  the  possibility 
lh.it  l.irge  CHOIS  m.iy  icsiill  when  ,i  lire.ikilown  ol  th.il  .ipproximalion  occurs. 


4.^ 


COMPARISONS  TO  EXPERIMENTAL  DATA 

Since  the  quantity  of  experimental  data  on  linewidth-vs-temperature  is  limited,  only 
a few  lines  will  be  considered  here.  It  is  important  to  recall  that  a line  shape  will  be  predicteil 
here  and  checked  with  the  experimental  results.  Only  normalization  to  room  temperature 
and  determination  of  the  low  temperature  intercept  are  necessary.  Using  the  tools  developed 
above,  consider  first  the  transition  in  Nd'*'^  YA(i  of  4p3/2(-)  ^ 41  j j ni  1 ) the  R2  ^ 1 

transition.  Referring  to  figure  5 showing  the  partial  energy  level  diagram  for  Nd'*'-^  \ A(i.  the 
following  is  found  for  the  general  expression,  before  simplification,  for  the  linewidth  of  the 
R2  Yi  transition: 


,R 


ZW(Rs- Y,)  = -IW'  +W 


R- 


(R)+  V w'^,.(Y)l  + 

7T  d21  ull 


1 = 2 


+ /Si^J  Y I ), 


R R 

whore  WJ:  and  VV’O  are  the  elastic  Raman  broadening  terms  for  the  Ra  and  \ i levels 

K T I 1 - I 

“ R /R 

respectively.  W^j j ( R ) is  the  downward  Raman  transition  from  R ->  to  R | . ) is  the 

upward  Raman  transition  from  ] to  Yj  ( i = 2 — .“i).  j(  R)  is  the  downward  direct-phonon 

transition  rate  from  R-i  to  R| . is  the  upward  direct-phonon  transition  rate  from 

Y]  to  is  any  multiphonon  transition  rate.  R-«)  is  the  inhomogeneous  random 

strain-induced  broadening  tor  the  Rs  level  and  ] ) is  the  strain  broadening  for  the  't'| 

level.  By  using  the  approximations  discussed  above,  combining  the  W^j  and  terms  into 
VV^  terms,  neglecting  the  direct  and  \\||^  terms,  and  neglecting  all  multiphonon  terms, 
the  following  is  used  as  the  expression  for  Ar't  Rs  ->■  't']  ): 


ArrtRs  - 'l  l ) = -i-  iW'U  +\v!J.  I +Ai^  (Rs)  + Ar^jY|). 

I TTKt  II  N-  ^l 


R R 

1 urihermore.  since  tlte  tunctional  dependence  with  temperature  of  W and  W. . is  identical 

K ■>  I I 

I independent  of  I ).  and  since  Ar’^(Rs)  and  Al^J.(^'| ) are  constants,  the  above  expressions  can 
be  further  reduced  to  the  extremely  simple  expression 


Ai^tRs  - V,  ) = jW'tRs  - Y|  ) + At^,^(Rs,  Y,  ). 


I sS) 


Now  all  that  is  reipiired  is  the  magnitude  of  \V*^(  R s 'i'  I ) ( related  to  the  constant  “r“  in 
equation  ( .sd ))  and  At',^1  R s.  | ):  rather  than  sevenfeen  difficult  to  measure,  coupled 
parameters  necessary  in  equation  t.s7).  fhe  number  of  unknowns  has  beeti  reduced  as  hn\  as 
possible  without  detailed  determinatioi.  , if  crystalline  wave  functions  and  matrix  elements. 
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Inn  W*^  = 0.  (59) 

I - 0 

the  meu!,ureinenl  of  R i " Y | ) at  low  temperatures  determines  the  residual  strain 
hroailenmg 

At't  R -)  “*  \ j ) = R T , Y I ).  ( 60 ) 

I - 0 

riieretbre,  h>  reciuirin^:  measurement  at  another  temperature  and  implicitly  assuming  that  Ar’j. 

IS  indepeiulent  ol  temperature  (obviously  violated  at  very  high  temperatures  causing  large  ionic 
displacements)  the  linewidth-vs-temperature  curve  for  the  R-i  -►  Y|  line  in  Nd'*’^  YAO  is  pre- 
dictetl,  I'or  convenience,  300°K  is  presumed  as  a reference  temperature  and  in  et]uation 
( 5S)  is  "normalized"  to  300°K.  F-'ormally,  a coefficient,  say  a,  is  introduced  and  fixed  by 
reij Hiring  that 

— W*^(300°K)=  1,  (61) 

7T 

or  if  the  ir  is  absorbed  into  the  tenn,  then  normalization  to  300°K  results  in 
( 7 , = T)  W'*^(300°K>, 

' IJ 

where  W ‘'(T)  is  given  by  b.quation  (54).  Therefore, 

Ai^(R2-*Y,)  = aW*^(T)  + Ai^j, 

with 

.1  = Ar'l  R 1 " 't  1 • “ ‘^'^s  ■ 

300°  K 

for  the  R-.-*Y|  curve,  the  low  temperature  interce]M  is  measureii  by  Kusliufr'*  to  lie 

0 X cni~*  and  the  300°K  linewidth  to  be  5.3  cm“^  fheretbre,  the  following  is  predieled  for 
the  liinetional  ilependence  with  temperature  of  the  R-«-*-Y  ] radiative  transition; 

R -1 Y I I = 4 5 cm” ' VV*^(  f)  + O.X  cm”  * (('3i 

When  I qii.ition  ((>3)  is  plotted  aiul  the  experimental  points  given  by  Kushida^”  are  includeit 
on  the  curve,  this  results  in  figure  14.  The  agreement  is  good  and  within  the  stated  experiniciit.il 
.iccuracy  of  ♦10'.'.  I lie  f.icl  tli.it  good  agreement  is  achieved  with  such  a simple  expression 
(fquation  (63)  was  plotted  with  a desk  calculator)  and  after  so  many  approximations,  is 
gr.ililyme  It  is  obseixed  tli.ii  the  tit  begins  (o  degrade  at  about  400°K  and  almost  10  errors 
.ire  seen  .it  470° K As  st.iied  earlier,  errors  .it  high  temperatures  are  not  unexpected  Since 
most  higher  onler  eltesis  .ire  deticiuleiit  on  the  partieul.ir  w.ixefunctions  involveil  (eg.  effects 

01  anh.irnionics  will  nalur.illy  be  wry  dependent  on  the  spati.il  form  of  ionic  wavefunctions), 

.1  prediction  of  when  Ingh-leinperature  errors  .ire  large  is  very  ilifficult.  It  will  be  observed 
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later  that  better  agreement  at  high  temperatures  is  achieved  with  some  levels  than  with 
others.  Also,  since  the  quality  of  the  fit  over  the  entire  curve  is  particularly  sensitive  to  the 
experimental  accuracy  of  the  measurement  of  linewidth  at  300°K  (and  somewhat  sensitive 
to  the  detennination  near  0°K ),  if  a large  error  occurs  experimentally  with  the  determina- 
tion of  the  linewidth  at  300°K,  the  shape  of  the  entire  curve  is  affected.  In  conclusion,  it 
can  be  said  that  the  agreement  between  the  theoretically  predicted  curve  and  the  experi- 
ment was  achieved  to  over  1 77°C  higher  than  most  Nd'*'^  YACi  lasers  operate). 

Consider  now  the  R | V'|  transition.  Proceeding  in  the  same  manner  as  used  for 

the  Ri  — Y|  transition,  the  low  temperature  intercept  is  measured  to  be  0.6  cm”'  and  the 
30()°K  point  is  4.0  cm”' . fherefore.  using  equation  (60). 

At'l R I Y[ ) = 3.4  cm”'  W*^(T)  + 0.6  cm”' . (64) 

This  result  is  plotted  along  with  experimental  points  in  figure  15  and  again  good  agreement  is 
realized.  The  maximum  error  occurring  at  445°K  is  nine  percent. 

Another  transition,  for  which  data  are  available,  is  the  primary  laser  line  in  Nd"'^'^  't'.AG 
This  line  is  composed  primarily  of  the  R-i  line  but  also  has  overlapped  from  the  R | 

(note  figure  1 5).  The  width  of  this  composite  line-vs-temperature  was  measured  by  Bell 
Laboratories'^  and  is  used  as  the  experimental  data  with  which  to  determine  the  (|uality  of 
fit  to  the  theory.  If  no.  consideration  is  given  to  the  fact  that  the  laser  line  is  actually  composed 
of  two  Lorentzian  lines,  but  the  two  measured  parameters  of  the  laser  are  used  with  eipiation 
(60),  a remarkably  good  fit  is  obtained.  Trom  this  data,  the  room  temperature  w kith  is  measured 
to  be  8.0  cm”'  and  the  low  temperature  width  is  0.5  cm”' . I'he  prediction  according  to 
et]uation  (58)  is  plotted,  along  with  the  data  point,  in  figure  16 

Aiz( laser  line)  = 7.5  cm”'  \v'^(T)  + 0.5  cm”'.  i(’5l 


I igiiie  14  The  ilicoretically  predicted  curve  lor  the  thermal  dependence  of  the  iinewidlh 
of  the  R-,  -*  Y|  radiative  transition,  showing  experimental  points  given  by  Kushida.'" 
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I igurc  1 5.  The  linewidth-vs-temperature  curve  predicted  for  the  R|  Y|  transition.'" 


F igure  16.  The  linewidth-vs-lenipcrature  curve  predicted  for  the  coiujxisite  lasei  line 


Here,  while  gmicl  aereenienl  is  observed,  u notieeahle  deviation  occurs  at  about  1 50°K  as  well 
as  the  expected  liigli-temperature  error.  The  deviation  around  1 50°K  is  attributed  to  tlie 
nature  ot  the  measurement  of  a "linew  idth"  composed  at  tvco  Lorent/ian  lines  ot  dilleren! 
magnitudes  separated  in  wavelength  added  to  a constant  strain-broadened  (iaussian  line. 
.\greernent  becomes  better  at  larger  temperatures  when  the  Lorent/ian  lines  become  much 
larger  than  the  C.aussian  lines. 

In  order  to  determine  the  transition  cross-section  and  branching  ratio  tor  the  laser 
line.  Kushida.  Marcos,  and  (ieuMc’"  decomposed  the  laser  line  into  components  of  Rt  V3 
and  R]  ^ '»'s  mote  figure  17)  ,\t  "room  temperature"  ( presumed  to  be  300°Kl  the  measured 
values  were  determined: 

ZWI  Rs  ^ '13)  I = .^.2  cm"’ 

* 300° K 

and 

At^lR]  ^ 't's)  =4.2  cm"’ 

300  K 

Notice  that  t his  correlates  extremeK  well  w ith  the  data  tor  the  R i ^ ')' I and  R [ * 't  | trans- 
ition where  linewidths  ot  .^.3  cm"’  and  4.0  cm  ’ were  observed  tthe  ditterence  ot  4.2  cm  ’ 
to  4.0  cm"’  is  percent  and  within  experimental  accuracv  ).  This  implies,  then,  that  the 
transitions  from  the  R levels  to  the  \ |.  ^ s and 't  ^ levels  are  broadened  almost  entireh  b\ 
the  R levels  alone,  and  the  contribution  to  the  broadening  due  to  the  'i  levels  via  an\'  mech- 
anism is  small.  This  deduced  fact  is  curious  though  unexplainable  w ithin  the  present  scope. 

.\n  attempt  to  explain  why  the  41  j | levels  are  not  significantly  thermally  broadened,  using 
the  present  tievelopment  w ill  not  be  attempted  here.  Therelore.  the  temperature  ilependence 
of  the  width  of  both  the  Rs  and  Rj  levels  has  been  explicitly  determined  and  can  be  used  in 
the  R ■ ransitions  to  the  /.  levels,  to  he  presented  next.  In  passing,  one  might  expect  a similar 
correlation  to  be  performed  on  the  residual  broadening^at  low  temperatures.  This  would  be 
meaningless,  though,  because  dillerent  samples  ol  Nil^''  'i  ,\(i  were  used  to  perlorm  the  ex- 
periment and  residual  strain  is  a function  of  sample  ciualiiy  and  can  exhibit  large  deviations 
sample  to  sample. 

The  next  radiative  transition  to  be  considereil  departs  from  the  R * nuiltiplets  and 
i.s  the  transition  from  the  TIN  i(  1 ) level  to  the  highest  4K)  s level,  the  Since  the  R| 
linewidth  is  already  determined,  all  that  is  needed  to  determine  the  linewidth  of  the  R | — /^ 
transition  is  the  temperature  depenilence  ot'the  /.^  level.  In  practice,  measurement  ot  the 
width  of  a single  level  is  not  practical  so  the  same  procedure  using  ei|uation  iPO)  w ill  be  used 
with  the  preceiling  data  Irom  just  Rj  used  to  establish  the  selt-cmisistency  ot  the  theorv.  Here, 
the  300  K linew  idth  is  found  to  be  ‘).2.x  cm"’  and  the  low -temperature  intercept  is  (1.0  cm"’ . 
This  results  in 

Ar’tRj  * 7.3  I = 3.2.^  cm”’  w’^(T)  + 6.0  cm"’ . (<>('1 

which  is  plotted  in  figure  lb  along  with  the  experimental  points.  The  agreement  is  generalK 
good  with  deviation  again  being  observed  at  high  temperatures:  10-percent  errors  are 
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observed  ;it  470°K.  Agiiin,  considering  tlie  approximations  made  to  simplitv  the  expressimi 
used  in  eejuation  (t)0),  and  taking  into  account  that  this  sliape  is  “predicted"  rather  than 
least-square  fit  to  tlie  data,  this  fit  is  considered  good.  If  the  magnitude  of  the  Raman  term 
in  the  R|  /.^  transition  cm"'  ) is  compared  with  the  previously  determined  value 

for  iust  the  R|  level.  3.4  cm"' . it  is  seen  that  here  the  width  of  the  R]  ^ Z5  transition  is 
composed  of  the  temperature-dependent  Raman-broailened  R|  level  and  the  temperature- 
independent  /.^  level.  .Actually,  there  is  a contribution  ot  the  strain  broadening  from  the 
R|  level  but  it  is  small  compared  to  b.O  cm"'  (from  equation  (<)4t  it  is  guaranteed  less  than 
0.6  cm"' ).  The  self-consistency  of  this  data  is  gratifying  (the  .^-percent  discrepancy  between 
the  two  values  is  within  the  experimental  accuracy ). 

In  the  preceding  section  it  was  stated  that  the  downward  direct-phonon  contribution 
to  broadening  was  combined  with  the  temperature-independent  strain  broadening  because 
for  most  cases  .vhen  l)'(to)  is  non-/,ero.  n(co)  + 1 1 and  is  independent  of  temperature. 

Since  it  is  really  unimportant  to  the  determination  of  the  linew  idth-vs-temperature  curve 
what  mechanism  causes  the  residual  broadening,  all  temperature-independent  terms  were 
combined.  I'or  the  level,  an  order  of  magnitude  increase  in  the  residual  broadening  is 
observed  anil  some  explanation  is  due.  Here,  the  Z;;  level  resides  .^37  cm  ' above  the  Z_j 
level  (figure  .‘s)  and  is  resonant  with  effective  phonons  contained  in  the  lattice  (within  n'(ecl). 
Since  n(.s37  cm"'  ) « 1 at  even  4.'^0°K.  this  results  in  an  estimate  of  what  the  coefficient  for 
direct  phonon  transitions  might  be  within  the  4lq  ->  nuiltiplet  (not.  that  is  out- 

siile  1)'(  CO  1) ; 

w'l^  (Z;;  Z_j)  s 6.0  cm"'  |n(.''3‘'  cm"'  kT)  + 1 1 l)'(.‘'3^  cm"'  ).  iO"i 

Hie  coefficient  of  6.0  cm"'  is  of  the  same  order  as  typical  coefficients  of  the  Raman  ex- 
pression namelv  ' 4 cm"' . This  direct-phonon  coefficient  is  ver\'  large  compared  to  the 
R ->  levels  already  examined  since  no  large  residual  broadenings  were  observed.  It  will 
be  shown  later  that  the  4l()  ->  levels  do  indeed  exhibit  extraordinariK’  large  direct-phonon 
transition  coefticients  whicli  result  in  a partial  breakdown  of  the  approximations  used  to 
derive  equation  (60>. 

The  final  curve  for  which  data  are  available  is  the  R|  * Z|  transition,  ll  the  low  tei^'p- 
erature  intercept  is  measured  to  be  0.5  cm"'  and  the  300''K  linewidth  is  measured  to  be 
I 3.2  cm"' . the  following  curve  is  derived  using  equation  (60): 

An(R,  - Z,  I = 11.'’  cm"'  \v'^(T)  -*■  0.5  cm"'. 

It  equation  (('bi  is  plotted,  and  data  points  entered . the  curve  m figure  l‘(  is  found.  I'he  tit 
is  bad,  .At  temperatures  ol  about  150  K,  errors  of  (lO  percent  are  observed  and  although  data 
points  .ire  taken  onK  to  320  K.  the  slopes  are  seen  to  be  radicalK  ditlerent.  This  poor  ft 
represents  .111  obiioiis  bre.ikdown  ol  cert.im  approximations  which  were  made.  It  will  turn 
out  th.it  the  sourie  ol  error  results  because  of  the  approximation  .is  to  the  shape  of  If'lo;), 
wherein  the  low  density  ot  si.iie  phonons  observed  from  ^5  cm"'  to  250  cm"'  were  ignored. 

I he  .ipproxim.ilion  w.is  nude  to  gre.itly  simplify  the  analy  tie  expression  lor  \\'*^  but  becomes 
iin.ihd  when  very  l.irge  direcl-phonon  transition  coelficients  occur  which  more  than  offset  the 
reduced  m.ignitude  ot  D’lcoi  tor  u;  • 250  cm"' . The  quantum  mech.mics  which  expl.nns  why 
the  4h,  s multiplets  h,i\e  .m  e\ti.iordm.irily  large  coefficient  is  not  understood  and  is  beyond 
the  present  scope  (eg.  group  theory  will  predict  only  qualitative  elfects;  qu.mtit.itive  results 
rei|uue  detei mm.il ion  ot  m.iiiix  elementsi. 
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Incorporation  of  tlic  direct  phonon  mechanism  will  now  be  pursued.  ,\s  stated  earlier, 
the  R|  level  broadenitis:  is  already  determined  and  in  tact,  since  the  residual  broadeninp  of  the 
R]  level  is  almost  identical  to  that  of  both  the  R]  and  /.  | levels,  all  that  is  to  be  dciermined 
is  the  temperature-dependent  level  broadening  of  the  Z|  level.  This  level  then  will  be  both 
Raman  and  upward,  direct  phonon  broadened.  If  the  Raman  broadening  of  the  Z|  level  is 
assumed  negligible,  the  upward  phonon  absorption  mechanism  remains,  and  the  following  can 
be  deduced ; 


/Si>[  R j ) = d.-'i  cm  ' \V*^(T ) + 0..S  cm  ' 


Ar^( Z 1 ) = X.2  cm  ' W (all  levels), 
I u 


wherein 


\\  (all  levels) 

' 300°  K 


350  400 


TEMPrRATUHE  ( K) 


I igiiio  I'l.  I lie  liiiew  kllli  As-teiiipoiatim.' ciiuc  picdiclod  tor  Ilic  R I • /.|  iransilioii. 


kitliciii  dioihuM  mccli.miMii 


H al!  multiplicative  coefficients  for  the  upward  transition  to  each  Stark  level  are  assumed  the 
same  ( there  is  no  reason  to  assume  this  except  for  convenience),  and  if  D’tto)  is  assumed  the 
same  for  all  upward  transitions  (this  is  rei|uired  because  the  detailed  shape  of  D'(tc)  fluctuates 
ftreatly),  the  folkiwinj;  can  be  written  for  the  linewidth  of  R|  / ] 

K I — Z I ) = d.5  cm“’  \V  ^(  f'l  + 8.2  cm~*  { | n(  1 24  cm“' : Tl  + 

+ n(  1 k7  em“^  ; T)  + n(3  1 1 cm~' ; T)  + n(848  cm“^ ; T)| 

Ini  134  cm”' ; 300‘’K)  + n(  107  cm"’ ; 300°K)  + 

+ 11(311  cm"' ;300°K) + 11(848  cm"' ;30()°K)1}.  '(P)i 

where  n(co;  T ) is  the  occupation  number  of  a phonon  of  fretiuency  to  at  temperature  T.  The 
plo*  t)f  etpiation  ((>81  ;ind  the  ilata  points  is  shown  in  fiiuire  20  and  excellent  agreement  is 
found.  Again,  the  partial  breakdown  of  the  theory  is  not  considered  critical  as  the  possible 
limits  to  the  validity  t)f  the  approximations  were  intiicateil,  and  also  the  tradeoff  of  simplicity, 
preilictivity , anil  consistency  is  significant. 


50  too  160  200  250  300  350  400  450 

TEMPERATURE  I Kl 


I igiiro  20  ihe  hitcwullh-v^-tcmpcrauiic  curve  lor  the  Rj  • /|  liaiiMlion  ulilamcil 
b\  emplot  mg  the  ditccl.  single  phonon  inechanisni '' 
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C ONC  Ll'SIONS 


As  stated  in  the  Introduetion.  tlie  objective  of  this  teehnieal  report  was  to  develop  a 
simple  compact  theory  to  predict  the  shape  ot  the  linewidth-vs-temperature  curves  lor  Nd^-’ 
VACi.  This  w'as  to  be  done  by  accomplishing  the  t'ollowing:  establish  the  internal  el'lestive 
mode  formalism  necessary  to  nKxiel  the  ion-lattice  interactions,  determine  a suitable  appros- 
imation  for  the  "effective"  phonon  spectra;  by  emplov  ing  this  approximation,  derive  compai  i. 
analytic,  single  multiplicative  parameter  expressions  for  the  broadening  mechanisms,  deter- 
mine the  ilominant  broadening  mechanism ; and  compare  the  predicted  results  w ith  the  experi- 
mental data.  The  previous  theories  had.  as  their  obiective.  the  verification  ot  the  vahditv 
ot  the  proposed  broadening  mechanisms  themselves  (direct,  Raman,  strain,  etc).  Ihe  succcns 
at  demonstrating  this  alone  has  been  excellent.  The  consir  imts  prohibiting  the  use  of  a 
Debye-based  theory  as  a means  of  predicting  the  thermal  variation  of  linewidlh  are  tour-told 
the  theoiy  of  Yen  employs  the  Debye  theoiy  for  phonons  which  is  not  valid  in  most  hosts 
(where  the  predominant  ion-lattice  interaction  is  via  optical  phonons)  and  rev|uires  a com- 
plicated curve-fit  procedure  to  determine  an  effective  ttjj  (ditferent  trom  the  resuliv  ot  bulk 
measurements);  the  functional  expression  for  the  Raman  contribution  is  nonanalvtic.  it 
requires  tlfe  measurement  of  many  points  on  many  linewidth-vs-temperature  curves  to  enable 
the  complicated  parametric  fit  to  be  accomplished;  and  finally,  and  most  signitlcant.  Ihe 
theory  of  Yen  does  not  "predict"  the  shape  of  the  curve. 

T he  internal  mode  theory  presented  here  alleviates  Ihe  above  constraints.  Here,  an 
extremely  simple  expression,  equation  (,sS)  is  u:wd  to  determine  the  thermal  variationv  ot 
Imevvidth  for  several  optical  transitions  between  multiplets  of  Nd'*’-''  'i  ACi  (for  which  data 
were  available).  The  R ^ Yj  and  R;  Yj  transitions  tit  well  with  the  predictions,  the 
expecletl  high  temperature  degradation  beginning  at  about  40(1  K.  This  again  in  attributed 
to  second-order  lattice  effects,  experimental  inaccuracies  m the  determination  ot  the  linew  idth 
at  ,^0()°K  and  near  0°K  used  to  determine  the  shape  of  the  curve,  and  increasing  Ihennal  oc- 
cupatiirn  of  low  energy  modes  resulting  in  small  errors  due  to  the  approximation  used  for 
D’(  co ) ( ie,  small  contributions  from  the  direct  mech.mism ).  ilhout  even  decomposition  of  the 
"laser  line."  the  fit  ot  this  transition  using  equation  (fiO)  was  excellent.  In  .iddilion  to  the 
expected  errors  at  high  temperatures,  Ihe  deviation  around  1 50  K was  attributed  to  errors 
resulting  Iron)  adding  two  Lorenl/ian  lines  ot  dillerent  m.igmiude  sep.irated  m wavelength 
added  to  a Ciaussian  line.  The  R]  ^ /.^  transition  also  tits  well  with  equation  ( 5S)  but  it 
was  clear  that  the  extraordinarily  large  residual  broadening  ot  ('  cm”*  is  due  to  direct,  spon- 
taneous emission  of  phonons,  since  Ihe  energy  g.ip  to  lower  levels  is  w ithin  D’lcc ).  I'his  is 
Ihe  first  indication  of  extraordinarily  large  ilirecl  phonon  coelTicienls  lor  nonradialive  tran- 
sitions within  the  4h)  s mulliplel.  When  the  same  proceilure.  using  e(|u.ilion  i(i0i.  w.is 
attempted  for  the  R|  • / | transition,  very  poor  results  were  ('btamed.  ll  is  Jc.ir  that,  in  this 
case,  the  ilirect  phonon  contribution  cannot  be  neglecleil.  When  mcliideil.  excellent  .igree- 
meiil  Is  again  achieved.  .As  slated  before.  Ihe  partial  breakilown  of  some  of  the  apprt'ximalions 
is  not  considereii  critical,  anil  indeed  they  were  anticipated. 

ll  can  be  stated  that  Ihe  objective  ot  this  theoretical  analysis  has  been  satisfied  ,i 
simple,  comp.icl  internal  moile  theory  predicting  the  shape  ot  Imewidlh-vs-temperature  curves 
lor  \d^'^  'i  \fl  has  been  developed. 
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